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ABSTRACT 



Context. The temperature and density structure of molecular cloud cores are the most important physical quantities that determine 
the course of the protostellar collapse and the properties of the stars they form. Nevertheless, density profiles often rely either on 
the simplifying assumption of isothermality or on observationally poorly constrained model temperature profiles. The instruments of 
the Herschel satellite provide us for the first time with both the spectral coverage and the spatial resolution that is needed to directly 
measure the dust temperature structure of nearby molecular cloud cores. 

Aims. With the aim of better constraining the initial physical conditions in molecular cloud cores at the onset of protostellar collapse, 
in particular of measuring their temperature structure, we initiated the Guaranteed Time Key Project (GTKP) "The Earliest Phases of 
Star Formation" (EPoS) with the Herschel satellite. This paper gives an overview of the low-mass sources in the EPoS project, the 
Herschel and complementary ground-based observations, our analysis method, and the initial results of the survey. 
Methods. We study the thermal dust emission of 12 previously well-characterized, isolated, nearby globules using FIR and submm 
continuum maps at up to eight wavelengths between 100//m and 1.2 mm. Our sample contains both globules with starless cores and 
embedded protostars at different early evolutionary stages. The dust emission maps are used to extract spatially resolved SEDs, which 
are then fit independently with modified blackbody curves to obtain line-of- sight-averaged dust temperature and column density maps. 
Results. We find that the thermal structure of all globules (mean mass 7 M Q ) is dominated by external heating from the interstellar 
radiation field and moderate shielding by thin extended halos. All globules have warm outer envelopes (14-20 K) and colder dense 
interiors (8-12 K) with column densities of a few 10 22 cm -2 . The protostars embedded in some of the globules raise the local tem- 
perature of the dense cores only within radii out to about 5000 AU, but do not significantly affect the overall thermal balance of the 
globules. Five out of the six starless cores in the sample are gravitationally bound and approximately thermally stabilized. The starless 
core in CB 244 is found to be supercritical and is speculated to be on the verge of collapse. For the first time, we can now also include 
externally heated starless cores in the L smm / L ho \ vs. r bo i diagram and find that T ho \ < 25 K seems to be a robust criterion to distinguish 
starless from protostellar cores, including those that only have an embedded very low-luminosity object. 

Key words. Stars: formation, low-mass, protostars - ISM: clouds, dust - Infrared: ISM - Submillimeter: ISM 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

** Partially based on observations carried out with the IRAM 30m 
Telescope, with the Atacama Pathfinder Experiment (APEX), and with 
the James Clerk Maxwell Telescope (JCMT). IRAM is supported 
by INSU/CNRS (France), MPG (Germany) and IGN (Spain). APEX 
is a collaboration between Max Planck Institut fiir Radioastronomie 
(MPIfR), Onsala Space Observatory (OSO), and the European Southern 
Observatory (ESO). The JCMT is operated by the Joint Astronomy 
Centre on behalf of the Particle Physics and Astronomy Research 



1. Introduction and scientific goals 

The formation of stars from diffuse interstellar matter is one of 
the most fundamental and fascinating transformation processes 
in the universe. Stars form through the gravitational collapse of 
the densest and coldest cores inside molecular clouds. The ini- 
tial temperature and density structure of such cores are the most 
important physical quantities that determine the course of the 
collapse and its stellar end product (e.g., |Larson|1969| |Penston 



Council of the United Kingdom, the Netherlands Association for 
Scientific Research, and the National Research Council of Canada. 



1 



R. Launhardt et al.: The structure of low-mass cloud cores 



19691 |Slml[T977l |Shu et al.|[T987| |Commercon eTaLl[20T0l ); 
however, deriving the physical properties of such cores from 
observations and thus constraining theoretical collapse models 
is a nontrivial problem. More than 98% of the core mass is 
locked up in H2 molecules and in He atoms that do not have 
a permanent electric dipole moment and that therefore do not 
radiate when cold. The excitation of ro-vibrational and elec- 
tronic states of H2 requires temperatures much higher than are 
present in these cores (e . g . , |B urton| 1 992| ) . This problem is usu- 
ally overcome by observing radiation from heavier asymmetric 
molecules, which are much less abundant, but have rotational 
transitions that are easily excited via collisions wit h hydrogen 
molecules (e.g., CO, CS; Bergin & Tafalla 2007). However, 
at the high densities (nu > 10°cm" J ) and low temperatures 
(T < 10 K; e.g., Craps reT"aL||2007| ) inside such cloud cores, 
most of these heavier molecules freeze out from the gas phase 
and settle on dust grains, where they remain 'invisible' (e.g., 



Ber gin & Lang er 1997; lCharnley|1997[|CalTelri et al.|1999[|Hn> 



Blant et al.||2010| ). The few remaining 'slow depleters' (e.g 
NH3 or N2H + , ftcrit(i-o) ~ 10 5 cm -3 ) or deuterated molecular 
ions that pr ofit from the gas p hase depletion of CO (e.g., N2D + 
and DCO + ; |Caselli et al.||2002| ) can be used to identify prestel- 
lar cores and infer their kinematical and chemical properties. 
However, the respective observational data, in particular for the 
deuterated species, are still sparse and their complex chemical 
evolution is not yet understood sufficiently well to derive reli- 
able density profiles. A very promising and robust alternative 
tracer for the matter in such cores is therefore the dust, which 
constitutes about 1 % of the total mass of the interstellar matter 
(ISM) in the solar neighborhood. 

Interstellar dust can b e traced by its eff ect on the attenua- 
tion of background light ( |Lada et aL]|1994 ), scattering of am- 
bient s tarlight ("cloudshine": [Foster & Go odman 2006; "core 



shine": |Steinacker et al.|2()l(), Pagani e t al.|2010| ), or by its ther- 
mal emission (e.g., Smith et al.|1979 ). In principle, it is prefer- 
able to use extinction measurements to trace the column density 
because they are to first order independent of the dust temper- 
ature (unlike emission). However, extinction measurements re- 
quire that at least some measurable amount of light from back- 
ground sources passes through the cloud, which is no longer 
the case at the highest column densities in the core centers 
(A^h > few x 10 22 cm -2 ). Another limitation of this technique 
is the inability to exactly measure the extinction toward individ- 
ual sources with with a priori unknown SEDs, and that it relies 
on statistically relevant ensembles of background stars, which 
restricts this method to fields with high background star density 
(i.e., close to the galactic plane) and limits the effective angu- 
lar resolution. Consequently, most attempts to derive the density 
structure of cloud cores from extinction maps have been done 
in the near-infrared (NIR) and have targeted less-opaque cloud 
cores and envelopes (A v < 20 mag; e.g.,|Alves et al.|2001[ [Lada 
|et al.|2004[ |Kandori et al.|2005| |Kainulainen et al.|2006| ). Going 
to wavelengths longer than optical or NIR, where the opacity 
is lower, also usually does not help much since the background 
stars also become dimmer at longer wavelengths. 

In addition to the conventional NIR extinction mapping tech- 
niques, mid and far-IR (MIR, FIR) shadows, or absorption fea- 
tures where the densest portions of cores are observed in absorp- 
tion against the background interstellar radiation field (ISRF), 
can be used to map structure at high resolution. Provided the 
absolute background level in the images can be accurately cal- 
ibrated, they provide good measurements of the line-of-sight 
(LoS) projected structure and column density in very dense re- 
gions. Such features have been used to study cores at 8, 24, and 



70 ^m ( [Bacmann efaLl[2000l |Stutz et al.|[2009a|b| |Tobin et al. 
|2010| ). One of the main results from these studies is that the 
LoS -projected geometry often drastically departs from spherical 
geometry both in the prestellar and protostellar phases. 

Observations of scattered interstellar radiation in the form 
of cloudsh ine, althou gh permittin g certain diagno stics of cloud 
structure ( |Padoan et al.||2006l |Juvela et al.| [2006 ), also do not 
trace the interiors of dense cores. Scattered MIR interstellar 
light can penetrate mo st low-mass cores producing coreshine 
( Steinacker et al. 2010). But, the effect requires the presence of 
larger dust grains for scattering to be efficient and is visible in 
only half of the cores (Paga ni et al.p OlO). Beside the advan- 
tage of an only weak dependence on temperature effects, scat- 
tered light is also sensitive to the 3D structure of the core and 
could allow to trace the actual density structure. But to use this 
information, the outer radiation field and the scattering phase 
function of the dust particles needs to be known precisely which 
poses a problem for many cores in the complex environment of 
star-formation regions. 

Hence, the thermal emission from dust grains remains as a 
superior tracer of matter in the coldest and densest molecular 
cloud cores where stars form. Consequently, most of our current 
information on the density structure of prestellar cloud cores and 
protostellar envelopes comes from submm and mm dust contin- 
uum maps (e.g.,|W ard-Thompso n et al.|1994|[T999}|Lau nhardt & 
lHenning|[T997l |Henning & Launhardt) 1 998 jlEvans et al.|2001[ 



Mot te & And re 2001; |Shirley et al.||2002[ |Andre et alT[[2 004; 
Kirk et al.|2005[|Kauffmann et al.|2008t|Launhardt et al.|2010p . 



Nevertheless, ultimately we will need to combine information 
from different tracers and methods to calibrate the dust opaci- 
ties as well as to derive quantities that cannot be derived from 
the dust emission alone (see discussion and outlook at the end of 
Sect.JT]). 

To relate the thermal dust emission at a given wavelength 
to the mass of the emitting matter, three main pieces of infor- 
mation are needed: the temperature of the dust grains, T d , their 
mass absorption coefficient, k v , and the gas-to-dust mass ratio. 
The mean opacity of the dust mixture is a function of grain op- 
tical properties, wavelength, and temperature (e.g., |He nning & 
Stognien ko||1996[ |Agladze et aLl|1996[ |Draine||2003HBoudet et 
|al.|2005| ) and may actually vary locally and in time as a function 
of the physical conditions. However, observational constraints 
of such opacity variations are very difficult to obtai n due to var- 
ious degeneracies and often insufficient data (e.g., |Shetty et al.| 
2009a|b| |Shirley et al.||20TTj |Kelly et aLpOH) . Therefore, the 



conversion of flux into mass is usually done by adopting one or 
the other 'standard' dust model that is assumed to be constant 
throughout the region of interest. 

The definition of a common temperature for the dust re- 
quires that the grains are in local thermal equilibrium (LTE), 
which is fortunately the case, to first order, in the well-shielded 
dense interiors of molecular cloud cores where collisions with 
H 2 molecules occur frequently and stochastic heating of indi- 
vidual grains (timescales r coo i « T h ea tX e -g-> by UV ph otons, 
does not play a significant role (e.g., Siebenmorg en et al.|1 992; 
|Pavly uchenk ov et al.|2012j ). However, with typical dust temper- 
atures inside such dense c ores being in the ra nge 8 - 15 K (e.g., 
di Francesco et al.||2007[ [Bergin & Tafalla 2007 ), the submm 



emission represents only the Rayleigh- Jeans tail of the Planck 
spectrum. Therefore, the dust temperature is only very poorly 
constrained by these data, which results in large uncertainties 
in the derived total masses and in the radial density profiles. 
E.g., at A = 1 mm, where the thermal emission of low-mass 
cores is practically always optically thin, more than twice the 
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amount of 8 K dust is needed to emit the same flux density as 
12 K dust. Hence, a dust temperature estimate of 10 ± 2K still 
leaves the mass and the (local) column density uncertain to a 
factor of two. The effect of unaccounted temperature gradients 
on the derived density profile is even more significant and leads 
to large uncertainties in the observational constraints on proto- 
stellar collapse models. Earlier attempts to reconstruct the dust 
temperature profiles of starless and star- forming cores from, e.g., 
SCUBA 450/850 /mi flux ratio maps or self -consistent radiative 
transfer modeling are discussed in Sect. 63 but remain observa- 
tionally poorly constrained and uncertain. 



This situation is currently improving as new observing ca- 
pabilities start to provide data that can potentially help to 
better constrain the dust temperatures. The Herschel Space 
Observatory with its imaging pho tometers that cove r the wave- 
lengths range from 70 to 500 fim ( [Pilbratt et al.|2010) provides 
for the first time the opportunity to fully sample the SED of the 
thermal dust emission from these cold object^] at a spatial res- 
olution that is both appropriate to resolve these cores and com- 
parable to that of the largest ground-based submm telescopes, 
which extend the wavelength coverage into the Rayleigh- Jeans 
regime of the SEDs. 



To make use of HerscheVs new capabilities and improve our 
knowledge on the initial conditions of star formation by over- 
coming some of the above-mentioned limitations, we have initi- 
ated the GTKP "The Earliest Phases of Star Formation" (EPoS). 
The main goal of our Herschel observations of low-mass pre- 
and protostellar cores in the framework of this project is to de- 
rive the spatial temperature and density structure of these cores 
with the aim of constraining protostellar collapse models. For 
this purpose, we selected 12 individual, isolated, nearby, and 
previously well-characterized cores and obtained spatially re- 
solved FIR dust emission maps at five wavelengths around the 
expected peak of the emission spectrum. This paper gives an 
overview of the observations, data reduction, analysis methods, 
and initial results on the temperature and column density struc- 
ture of the entire low-mass sample. Initial results on one of the 
sources (CB 244) were already published by |Stutz et al. ( 2010] ). 
A first detailed follow-up s tudy of another source from this sam- 
ple (B 68) is published by |Nielbock et al.| ( [20T2] ). Initial results 
on the high- mass cores observed i n the frame work of EPoS ar e 
publ ished by|Beuther et al.| ( [20T0|[20T2] ), |Henning et al.| ( [20T0] ), 
and |Linz et al.| ( 2010|). An o verview on the high-mass part of 
EPoS is given in Ragan et al.](|2012|). 



This paper is organized as follows. Section|2] describes 
the target selection criteria and the sample of target sources. 
Section]!] describes the Herschel and complementary ground- 
based observations and the corresponding data reduction. In 
Sect. [4] we introduce our method of deriving dust temperature 
and column density maps from the data. Section]?] gives an 
overview of the initial results from the survey, and in Sect.|6] 
we discuss the uncertainties and limitations of our approach and 
compare our results to earlier work by other authors. Section]?] 
summarizes the paper. 
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Fig. 1. All-sky map (gnomonic projection), showing the mean 
stellar K-band flux density distribution of the Milky Way 
(grayscale, derived from COBE-DIRBE NIR all-sky maps 2 ) and 
the location of our target sources (see Table]!]). 
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Fig. 2. IRAS lOOyimi dust continuum emission map, showing 
the p Oph region and the loaction of two of our target sources 
(B 68 and CB 68). The plane of the Milky Way becomes visible 
at the lower left corner of the map. Compare to Figs.[T8]and|A.7 



2. Sources 



Based on the results of earlier studies (e.g., Launhardt & 



Henning 1997 Launhardt et al. 1997 Henning & Launhardt 



1998[ |Stutz et al.||2009aj |Launhardt et al.||2010| ), we have se 
lected 12 well-isolated low-mass pre- and protostellar molec- 
ular cloud cores in regions of exceptionally low cirrus confu- 
sion noise. In order to obtain deep maps at 100//m, which is 
essential for a precise estimate of the dust temperature, the ab- 
solute background levels and the point source confusion noise 
(PSCN) were important selection criteria. For this purpose, con- 
tributions from spatial flu ctuations of the extr agalactic back- 
ground were derived from Negrello et al. (2004). Estimates of 
the Galactic cirrus confusion noise are based on the ISOPHOT 
confusion noise measurements, scaled down in the power spec- 



1 The peak of the flux density of a blackbody with T = 10 - 20 K 
is at A « 290 - 145 iim. The peak of the energy density (vS v ) of the 
optically thin globules at these temperatures is at only slightly shorter 
wavelengths (230 - 130//m). 

2 http://lambda.gsfc.nasa.gov/product/cobe/dirbe_overview.cfm 
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Table 1. Source list 



Source 


Other 


R.A., Dec. (J2000) 


Region 


Dist. 


Ref. 


Evol. 


Ref. 




names 


[h:m:s, °:':"] 




[pc] 




class 




CB4 




00:39:03, +52:51:30 


Cas A, GB 


350 + 150 


18 


starless 


1 


CB 6 


LBN 613 


00:49:29, +50:44:36 


Cas A, GB 


350 + 150 


18 


CI. I 


2 


CB 17 


L 1389 


04:04:38, +56:56:12 


Perseus, GB 


250 + 50 


2,18 


starless^ 


2 


CB 26 


L1439 


05:00:09, +52:04:54 


Taurus -Auriga 


140 + 20 


2,11 


starless 


2,3,4 


Cr$ Z / 


L, IjIZ 


Uj.U4.Uy, +JZ.43.Uo 


Taurus - Auriga 


1 /in _i_ on 
14U + ZU 


111 


starless 


1 A 

1,4 


BHR12 


CG30 


08:09:33, -36:05:00 


Gum nebula 


400 + 50 


2,10 


CI. C 


2,10 


CB 68 


L146 


16:57:16,-16:09:18 


Ophiuchus 


120 + 20 


2,11 


C1.0 


2 


B68 


L57, CB 82 


17:22:35, -23:49:30 


Ophiuchus, Pipe nebula 


135 + 15 


15,16,17 


starless 


5 


CB130 


L507 


18:16:15, -02:32:48 


Aquila rift, GB 


250 + 50 


12,14 


CI. d 


2 


B335 


CB199 


19:36:55, +07:34:24 


Aquila, isolated 


100 + 20 


8 


C1.0 


2,7 


CB 230 


L1177 


21:17:39, +68:17:32 


Cepheus flare, GB 


400 + 100 


2,13 


CI. I 


2 


CB 244 


L1262 


23:25:47, +74:17:36 


Cepheus flare, GB 


200 + 30 


2,13 


CI. o e 


2,9 



Notes. (a) CB 17: additional low-luminosity Class I YSO 25" from starless core. {b) CB 26: additional Class I YSO 3.6' south-west of the starless 
core. (c) BHR 12: additional Class I core ~20" south of Class source. {d) CB 130: additional low-mass prestellar core ~30" west and Class I YSO 
~15" east of Class core. (e) CB 244: additional starless core ~90" west of Class source. 

References. (1) this paper; (2)|Launhardt et aLlpOTO]); (3)[Launhardt & S argent) pom]); (4)|Stutz et al.]p009a|); (5)|Alves et al.lpOOT]); (7)|Stutz 
leTaLl([2008l); (8 ) |Qlofsson & Qlofsson| <|2009l>; (9) |Stutz et al. 142010^ (10) Bour ke et al.| ( |1995|); (ll)|Loinard et al.|<2011); (12)|Launhardt & 
Henni ngl([T997|>; (13)|Kun|(|1998| ); (14) |Straizys et al.| ( |2003| ); (15) |de Geus et al.| ( |1989| >; (16) |Lombardi et al.| ( |2006| )r(17) |Alves & Franco N2007] ); 
(18) |Perrot & Grenier| ( |2003] r~ 



trum to the resolution of Herschel/PACS ( [Kiss e t al. 2005}. Our 
selected globule sources have typical lOOyurn background levels 
of lmJy/ D// or lower and PSCN< 1 mJy/beam. For compari- 
son, typical lOOyum background levels within the Taurus cloud 
are about 2 - 4mJy/ D// , within Ophiuchus about 5mJy/ D// , and 
reach values of > 100mJy/ D// in regions of high-mass star for- 
mation and infrared dark clouds (IRDCs). This specific source 
selection strategy enabled us to obtain deeper 100 yum maps than 
most large-area surveys and to derive robust dust temperature 
estimates. 

Our target list contains only established and previously well- 
characterized sources and does not cover regions with unknown 
source content. All sources are nearby (100- 400 pc, with a 
mean distance of 240 ± 100 pc), have angular diameters of 3' 
to 6', linear sizes between 0.2 and 1.0 pc, and total gas masses of 
1 to 25 M (see Table[6]). Coordinates, distances, spectral (evo- 
lutionary) classes, and references are listed in Table[T] FigurefT] 
shows the distribution of the target sources on the sky. Figure|2j 
illustrates that even those globules in our source list that are 
loosely associated with larger dark cloud complexes (like p Oph) 
are still isolated and located outside the regions of high extinc- 
tion and source confusion. 

Seven out of the 12 target globules were known to con- 
tain starless cores ([Ward-Thomps on et al.|1 994). Of these, only 
CB 17 was already shown before to be prestellai^] in nature 
(Pavlyuchenkov et al. 2006 ), while the star-forming potential of 
the other sources is not known. CB 17 contains, in addition to 
the prestellar core, a low-luminosity Class I young st ellar object 
(YSO) at a projected separation of -25" (-6000 AU; Launhardt 
et al. 2010). CB 26 contains, in addition to the starless core ( [Stutz 
et al.|2009a| , a well-st udied Class I YSO at a projected separa- 
tion of ~3.6' (~0.15pc;|Lau nhardt & Sargent 2001 ; Stecklum et 
|al. 120041 |Launhardt et al.|2009[|Sauter et al.|2009| . Five out of 
the 12 target globules were known to host Class protostars. Of 
these, three cores contain additional sources with other evolu- 
tionary classifications. CB 130 contains an additional low-mass 
starless core as well as a Class I YSO. BHR 12 contains an addi- 



tional Class I core at a projected separation of -20" (-8000 AU). 
CB 244 contains an addition al starless core at a projected sepa- 
ration of -90" (~18000AU; |Launhardt et al.|[20i0l|Stutz et al.| 
|20T0| ). Two out of the 12 target globules (CB 6 and CB 230) are 
dominated by embedded Class I YSOs, i.e., the extended emis- 
sion in them supposedly arises from remnant post-collapse en- 
velopes. Of these, CB 230 is known to be a binary source with a 
projected separation of -10" (-4000 AU; |Launhardt et al.|2010l ). 



We have also re-evaluated the distance estimates toward all 
globul es, confirming the ea rlier used distances for 7 sources 
(e.g., Launhardt et al. 2010), and adjusting the distances for 5 
sources. CB 4 and CB 6, which we earlier associated with the 
so-called "-12km/s" HI clouds at 600- 800 pc, are unlikely 
to be that far away, as suggested by the lack of foreground 
stars within a 3' diameter area toward the cores. Despite their 

vlsr 12 km/s, they are more likely a ssociated with the Cas A 

dark clouds in Gould's Belt at ~ 350 pc ( [Perrot & Grenier|2003| ). 
CB 68 is associated with the p Oph dark clouds, for which we 
adopt the new precis e trigonometric distance of 120pc from 
|Loinard et al.| ( [201 1| ). B 68 is located within the Pipe nebula 
and somewhat farther away from p Oph. Various d irect and indi- 
rect estimates suggest a distance of 135 pc ±15 pc (de Geus et al. 



198^ [Lombardi et al.|2006[|Alves & Franco|2007| ). For CB 130, 
which is associated with the Aquila R ift clouds, we adop t the 
distance estimate of 250 pc ±50 pc from Straiz ys et al.| ([2003 ). 



3 We use the term "prestellar" only for those starless cores that have 
been shown to be gravitationally unstable or at the verge of collapse. 



3. Observations and data reduction 

The Herschel FIR continuum data are complemented by submm 
dust continuum emission maps at 450 yum, 850 yum, and 1.2 mm, 
obtained with ground-based telescopes. Observations and reduc- 
tion of both the Herschel and complementary ground-based data 
are described in the following two subsections. Representative 
general observing parameters are listed in Table[2| and the ob- 
servations of individual sources are summarized in Tabled 
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Table 3. Summary of FIR and submm observations 



Source 


PACS 


SPIRE 


450yL/m 


850yL/m 


1.2 mm 




ObsID 


ObsID 








CB4 


134221614(5,6) 


1342188689 


- 


- 


Mil 


CB6 


134222272(3,4) 


1342189703 


L10 


L10 


L10/M11 


CB17 


134219100(8,9) 


1342190663 


- 


L10 


L10/M11 


CB 26 


134221739(7,8) 


1342191184 


L10 a 


L10 a 


L10/M11 


CB 27 


134219197(5,6) 


1342191182 


DF08 


DF08 


Mil 


BHR12 


13421985(39,40) 


1342193794 


L10 


L10 


LIO 


CB 68 


134220428(6,7) 


1342192060 


L10 


L10 


LIO 


B68 


134219305(5,6) 


1342191191 


DF08 


DF08/N12 


B03 


CB130 


134220598(2,3) 


1342191187 


L10 


L10 


LIO 


B335 


134219603(0,1) 


1342192685 


L10 


LIO 


Mil 


CB 230 


134218608(3,4) 


1342201447 


L10 


LIO 


L10/K08 


CB 244 


134218869(4,5) 


1342199366 


L10 


LIO 


LIO 



Notes. {a) The SCUBA maps of CB 26 cover only the western core with the embedded Class I YSO, but not the eastern starless core. 

References. (L10) |Launhardt et al.|f2010| >; (Mil) dedicated MAMBO- 2 observations, described in this pap er (see Sect.[3^2|) ; (DF0 8) |Di Francesco| 
|eTaI] ( [2008l SCUB A Legacy Catalogues); (B 03) |Bianchi et al.| < [20Q3] SEST-SIMBA observations); (N12) |Nielbock et al'lpOHj APEX-Laboca 
observations); (K08) |Kauffmann et al.| f2008 IRAM-MAMBO-2 observations) 



Table 2. General observing parameters 



^0 


Instr. 


(HPBW> 


(map size) 


(rms) 


[jim] 




[arcsec] 


[arcmin] 


[mJy/beam] 


100 


PACS 


l.V 


10 


3 


160 


PACS 


\\2 a 


10 


13 


250 


SPIRE 


18.2 a 


16-20 


15 


350 


SPIRE 


25.0 a 


16-20 


17 


450 


SCUBA 


8.5-9.5^ 


2 


200 


500 


SPIRE 


36.4 a 


16-20 


14 


850 


SCUBA 


14.4-15.1* 


3 


20 


870 


LABOCA 


19.2 C 


40 


8 


1200 


MAMBO-2 


11 


5-10 


6 



Aniano et al.|j201lj > (Z?) |Launhardt et al.| ( |2010| ) (c) |Nielbock 



Notes. (a \ 
|eTaLl ( [20T2F 



3. 1. Herschel FIR Observations 

FIR continuum maps at five wavelength bands were obtained 
with two different instruments on board the Herschel Space ob- 
servatory: PACS at 100 and 160 /mi, and SPIRE at 250, 350, and 
500 jim. 

3.1.1. PACS data 

All 12 sources were observed wi th the Photodetector Array 
Camera and Spectrometer (PACS; Poglitsc het al.|2 010) in the 
scan map mode at 100 and 160 /mi simultaneously. For each 
source, we obtained two scan directions oriented perpendicular 
to each other to eliminate striping in the final combined maps of 
effective size ~ 10' x 10'. The scan speed was set to 20" sec -1 , 
and a total of 30 repetitions were obtained in each scan direc- 
tion, resulting in a total integration time of ~ 2.6 firs per map. 
We emphasize that the accurate recovery of extended emission 
is the main driver for the following discussion and exploration of 
reduction schemes. The PACS data at 100 fim and 160 fim were 
processed in an identical fashion. They were reduced to level 1 
using HIPE v. 6.0.1196 ( [Qtt|[20T0| ), except for CB4 (processed 
with HIPE v. 6.0.2044), CB26 (processed with v. 6.0.2055), and 
CB6 (processed with v. 7.0.1931). The rationale for using dif- 
ferent versions of HIPE for different sources is that the latter 
three sources were observed much later (OD 660-770) than the 



other sources (OD 230-51 1) and the data were reduced with the 
respective latest HIPE version and calibration tree that was avail- 
able at that time. Since processing of the PACS scanning data is 
very time-consuming, a re-reduction of all data with every new 
HIPE version with modifications relevant for our data is not fea- 
sible. Instead, we re-reduced the data for one source with the lat- 
est version of HIPE to verify how much the data reduction and 
post-processing modifications affect our results on the derived 
dust temperature and column density maps. Since this compari- 
son, which is presented and discussed in Sect. 6.1 showed that 
the effects are well within the range of the other uncertainties, 
we do not re-reduce all data, but incorporate these effects into 



our uncertainty assessment in Sect. 6.1 

Apart from the standard reduction steps, we applied '2 nd 
level deglitching' to remove outliers in the time series data 
{'time-ordered' option) by applying a clipping algorithm (based 
on the median absolute deviation) to all flux measurements in 
the data stream that will finally contribute to the respective map 
pixel. For our data sets we applied an 'nsigma ' value of 20. 

After producing level 1 data, we have generated final level 2 
maps using two different methods: high-pass filtering with phot- 
project (within HIPE) and Scanamorphos (Roussel [2012| ), which 
does not use straightforward high-pass filtering, but instead ap- 
plies its own heuristic algorithms to remove artifacts caused by 
detector flickering noise as well as spurious bolometer temper- 
ature drifts. Based on the high-pass median- window subtrac- 
tion method, the photproject images turned out to suffer from 
more missing flux and striping than the Scanamorphos images. 
Because the correct recovery of extended emission is critical for 
accurate temperature mapping, high-pass filtering was clearly 
disadvantageous for our science goals. We also note that in a 
previous reduction, we found that the MADmap^] processing of 
our brighter sources (e.g., B 335) used to produce large artifact^] 
in the final maps due to the relatively bright central protostar; 
therefore we do not utilize this reduction scheme in this work. 

For these reasons we have decided to use Scanamorphos v. 9 
as our standard level 2 data reduction algorithm. A comparison 



4 http://herschel.esac.esa.int/twiki/pub/Pacs/PacsMADmap/- 
UM_MADmapDOC.txt 

5 This problem has eventually been solved in the most recent 
MADmap implementation. 
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of final maps processed with different Scanamorphos v. 9 options 
showed that the "galactic" option recovered the highest levels 
of extended emission and was thus the best-suited for our data 
set and scientific goals. We use a uniform final map pixel scale 
of 176 pix" 1 at 100 /mi and 372 pix" 1 at 160 yum for all ob- 
jects. Furthermore, these data were processed including the non- 
zero-acceleration telescope turn-around data, with no additional 
deglitching ('no glitch' setting). 

3.1.2. SPIRE data 

All 12 sources were also observed at 250, 350, and 500 yum with 
the Spectral and Photometric Imaging Receiver (SPIRE; Griffin 
|et al.|2010| ). The scan maps of each source were obtained simul- 
taneously at all three bands at the nominal scan speed of 30"/sec 
and over a (16' - 20' ) 2 area, resulting in a total integration time 
of 10 - 16 mins per map. These data were processed up to level 1 
with HIPE v. 5.0.1892 and calibration tree 5.1, which was the 
most u p-to- date version at the time all data became available. 
In Sect. |6.1| we show that using the newest HIPE v. 9.1.0 only 
leads to nonsignificant minor changes in the SPIRE maps that 
are well below all other uncertainties. For this reason, we consid- 
ered it unnecessary to re-reduce the data. Up to level 1 (i.e., the 
level where the pointed photometer timelines are derived), the 
steps of the official pipeline (POF5 -pipeline. py, dated 2.3.2010) 
provided by the SPIRE Instrument Control Center (ICC) have 
been performed. The level 1 frames were then processed with 
the Scanamorphos software version 9 ( |Roussel|2012| ). All maps 
were reduced using the 'galactic' option. The final map pixel 
sizes are 6, 10, and 14" at 250, 350, and 500 yum, respectively. 

Photometric color corrections for both PACS and SPIRE 
data, that account for the difference in spectral slopes between 
flux calibration (flat spectrum) and actually observed source 
spectrum, are applied only in the data analysis and are described 
inSect.fO 



3.2. (Sub)mm continuum Observations 

Complementary submm dust continuum emission maps at 
450 yum, 850 yum, and 1.2 mm from ground-based telescopes are 
available or were obtained through dedicated observations for all 
12 sources. References to the existing and published data used 
here are given in Table[3] 

Additional dedicated maps of the 1.2 mm dust emission 
from 6 sources were obtained with the 117-pixel MAMBO-2 
bolometer array of the Max-Planck-Institut fur Radioastronomie 
( [Kreysa et al.||1999] ) at the IRAM 30 m-telescope on Pico 
Veleta (Spain) during two pool observing runs in October and 
November 2010. The mean frequency (assuming a flat- spectrum 
source) was 250 GHz (Aq ~ 1.2 mm) with a half power (HP) 
bandwidth of -80 GHz. The HP beam width (HPBW) on sky 
was 11 ", and the FoV of the array 4' (20" pixel spacing). 
Weather conditions were good, with zenith optical depths be- 
tween 0.1 and 0.35 for most of the time and low sky noise. 
Pointing, focus, and zenith optical depth (by sky dip) were mea- 
sured before and after each map. The pointing stability was bet- 
ter than 3" rms. Absolute flux calibration was obtained by ob- 
serving Uranus several times during the observing runs and mon- 
itored by regularly observing several secondary calibrators. The 
flux calibration uncertainty is dominated by the uncertainty in 
the knowledge of the planet fluxes and is estimated to be -20%. 
The sources were observed with the standard on-the-fly dual- 
beam technique, with the telescope scanning in azimuth direc- 



tion at a speed of 8"/s while chopping with the secondary mirror 
along the scan direction at a rate of 2 Hz. Most sources were 
mapped twice, with different projected scanning directions (ide- 
ally orthogonal) and different chopper throws (46" and 60") to 
minimize scanning artifacts in the final maps. For one source 
(CB 27), we could obtain only one coverage. Effective map sizes 
(central pixel coverage) were adapted to the individual source 
sizes and were in the range 5' - 10', resulting in total mapping 
times per coverage of 40 min to 1.5 hrs. 

The raw data were reduced using the standard pipeline 
"mapCSF" provided with the MOPSIC softwar^] Basic reduc- 
tion steps include de-spiking from cosmic ray artifacts, cor- 
related signal filtering to reduce the skynoise, baseline fitting 
and subtraction, restoring single-beam maps from the dual-beam 
scan data using a modified EKH algorithm ( [Emerson et al.| 
|1979| ), and transformation and averaging of the individual hori- 
zontal maps into the equatorial system. Correlated skynoise cor- 
rection was performed iteratively, starting with no source model 
in the first run. The resulting map is smoothed and the region 
with source flux is selected manually to represent the first source 
model, which is then iteratively improved in 20 successive corre- 
lated skynoise subtraction runs. Since our sources are faint, we 
follow the recommendation to apply this algorithm only up to 
level 1, except for the brightest source B335, where we go up to 
level 3. A third order polynomial baseline fit to the time-ordered 
data stream as well as first order fits to the individual scan legs 
are subtracted from the data, after masking out the region with 
source flux, to correct for slow atmospheric and instrumental 
drifts. Despite some unsolved electronic problems during these 
observing runs, which limited the efficiency of the correlated 
skynoise subtraction, the resulting noise levels in the final maps 
are ~6mJy/beam, except for CB 27, where it is ~11 mJy/beam 
because only one coverage could be obtained. 

The assumption of a flat-spectrum source in the definition 
of the nominal wavelength of the MAMBO detectors in com- 
bination with the broad bandpass made it necessary to ap- 
ply a "color correction". In contrast to the Herschel FIR ob- 
servations, (sub)mm wavelengths are sufficiently close to the 
Rayleigh- Jeans regime and the emission is optically thin for all 
our sources, such that the slope of the SED follows S v oc A~ 4 
(for the dust opacity submm spectral index = 2). Therefore, 
we follow Schnee et al. (2010) and correct the reference wave- 
length for all MAMBO data to 11 00 /mi. Color corrections for 
the SCUBA data are negligibly small (cf. [Schnee et aL]|2010| ) 
and were not applied. 

We also verified the pointing in the final maps by comparing 
the peak positions with interferometric positions where possi- 
ble (CB 26, CB 68, CB 230, CB 244, B 335, BHR 12) and found 
that deviations are in all cases smaller than 2", hence confirming 
the pointing stability and making additional pointing corrections 
unnecessary. The final maps of those sources for which we had 
both old and new 1 .2 mm data (CB 6, CB 17, CB 26, and CB 230) 
were generated by weighted averaging, after adapting beam sizes 
and correcting small pointing offsets. 

Additional complementary NIR extinction maps as well as 
different molecular spectral line maps were also obtained with 
the aim of studying the relation between density, temperature, 
and dust properties on the one side, and gas phase abundances 
and chemistry on the other side. These latter data will be de- 
scribed and analyzed in forthcoming papers and are not pre- 
sented here. 



6 Created and continuously updated by R. Zylka, IRAM, Grenoble; 
see http://www.iram.es/IRAMES/mainWiki/CookbookMopsic 
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Table 4. DC-flux levels in the Herschel data 



Source 


R.A.,Dec. fl 


fDC (O-f) 


fDC (Cf) 


fDC (Cf) 


fDC (O-f) 


fDC (O-f) 




(J2000) 


lOOyL/m 


160 /mi 


250 //m 


350yL/m 


500 //m 

^Jy/ D,/ ] 




[h:m:s, °:':"] 


^Jy/ D,/ ] 


^Jy/ D "] 


[//Jy/ D "] 


^Jy/ D,/ ] 



CB4 
CB 6 
CB17 
CB 26 
CB 27 
BHR12 
CB 68 
B68 
CB130 
B335 
CB 230 
CB 244 



00:38:55, 
00:49:02, 
04:04:14, 
05:00:28, 
05:04:16, 
08:09:10, 
16:57:07, 
17:22:36, 
18:16:30, 
19:36:48, 
21:17:00, 
23:24:55, 



+52:56:06 
+50:42:34 
+56:53:07 
+52:00:36 
+32:38:19 
-36:07:44 
-16:14:27 
-23:44:11 
-02:29:23 
+07:30:37 
+68:14:35 
+74:23:01 



8.4E2 
1.1E3 
1.3E3 
7.8E2 
1.0E3 
1.3E3 
1.3E3 
9.1E2 
1.4E3 
1.5E3 
1.3E3 
1.3E3 



(4.4E1) 
(3.0E1) 
(4.2E1) 
(4.7E1) 
(4.7E1) 
(5.7E1) 
(5.4E1) 
(8.4E1) 
(7.1E1) 
(4.3E1) 
(4.8E1) 
(4.0E1) 



5.9E2 
4.6E2 
7.5E2 
7.7E2 
6.7E2 
6.7E2 
5.1E2 
8.9E2 
8.0E2 
7.2E2 
5.1E2 
5.1E2 



(1.0E2) 
(3.7E1) 
(4.5E1) 
(1.1E2) 
(1.3E2) 
(8.5E1) 
(6.9E1) 
(7.2E1) 
(1.5E2) 
(7.9E1) 
(1.8E2) 
(7.5E1) 



2.8E2 
1.8E2 
2.6E2 
1.4E2 
1.9E2 
3.1E2 
2.2E2 
1.2E2 
3.0E2 
1.5E2 
1.8E2 
3.0E2 



(3.9E1) 
(2.2E1) 
(3.1E1) 
(2.7E1) 
(2.8E1) 
(2.9E1) 
(3.5E1) 
(2.9E1) 
(7.4E1) 
(2.1E1) 
(4.3E1) 
(2.7E1) 



1.4E2 
8.7E1 
1.3E2 
6.8E1 
9.4E1 
L7E2 
1.3E2 
5.9E1 
1.6E2 
8.7E1 
1.4E2 
1.7E2 



(2.1E1) 
(1.1E1) 
(1.6E1) 
(1.0E1) 
(1.7E1) 
(L3E1) 
(2.0E1) 
(L4E1) 
(3.9E1) 
(9.0E0) 
(7.1E1) 
(1.3E1) 



6.9E1 
3.8E1 
6.5E1 
3.3E1 
5.0E1 
7.7E1 
4.2E1 
1.3E1 
7.7E1 
3.9E1 
7.4E1 
7.7E1 



(1.0E1) 
(5.0E0) 
(7.0E0) 
(7.0E0) 
(7.0E0) 
(5.0E0) 
(1.1E2) 
(5.0E0) 
(1.7E2) 
(6.0E0) 
(2.8E1) 
(7.0E0) 



Notes. (a) Center of the DC reference field (see Sect.Ej}. 



4. Modeling approach 

4. 1. Strategy and data preparation 

The calibrated dust emission maps at the various wavelengths 
were used to compile spatially and spectrally resolved data cubes 
that cover the full extent of these relatively isolated sources on 
both sides of the peak of their thermal SEDs. Ideally, one would 
want to compare such data to synthetic maps from radiative 
transfer models convolved with the respective beams (forward- 
modeling, virtual observations) in order to constrain the phys- 
ical properties of the sources. However, to avoid circular aver- 
aging with its well-known caveats, one would need 3-D model- 
ing to account for the complex structure present even in these 
relatively simple and isolated sources. Furthermore, not all of 
the observed featu res ( like , e.g ., the core-envelope temperature 
contrast, see Sects. ^6 and |5.7| ) may be easily reproducible with 
existing self-consistent models, making the fully self-consistent 
forward-modeling approach very time-consuming. This will be 
dealt with in forthcoming papers modeling individual cores. 

For this survey overview paper, we therefore take a simpler 
and more direct approach, giving up some of the spatial infor- 
mation in the short-wavelength maps and directly recovering 
and modeling beam and LoS optical-depth-averaged (hereafter 
LOS-averaged for short) SEDs, thus introducing as few as pos- 
sible model-dependent assumptions into our analysis. For this 
purpose, the calibrated maps are prepared in the following way: 

1 . all maps are registered to a common coordinate system, 

2. pointing corrections are applied where necessary (see be- 
low), 

3. all maps are converted to the same physical surface bright- 
ness units (here Jy/ D// ) and extended emission calibration 
corrections are applied where necessary (see below), 

4. background levels are determined and subtracted from the 
maps (see below), and finally, 

5. all maps are convolved to the SPIRE 500 yum beam (FWHM 
3674). 

PACS pointing corrections: Since the Herschel spacecraft abso- 
lute 1 o~ pointing error was ~2" and the PACS and SPIRE maps 
were not obtained simultaneously, random relative pointing er- 
rors between PACS and SPIRE maps can add up to ~3"-4", 



which is abo ut half the PACS 100 yum FWHM beam size ( [Aniano 
et al.|20lT ). This could severely hamper the usefulness of flux 
ratio (SED) maps, in particular around compact sources, if not 



corrected properly. For the long-wavelength SPIRE maps with 
beam sizes ~18" -36" and no readily available astrometric ref- 
erence, the pointing errors are fortunately negligible. For the 
PACS maps we utilize the fact that many lOO/mi point sources 
are also detected in the Spitzer MIPS 24 jim images. We select 
in each field three stars detected in both the PACS lOO/mi and 
MIPS 24 /mi images to align the PACS astrometry in both the 
lOOyurn and 160yum images to the Spitzer images. The pointing 
corrections were found to be < 3" in all cases, i.e., small com- 
pared to the SPIRE 500 yum beam. 

Surface brightness calibration corrections: Since the PACS im- 
ages are calibrated to Jypix -1 , the conversion to Jy/ D// does not 
involve an assumption of the beam sizes. The SPIRE data, on 
the other hand, are calibrated to Jy beam -1 ; to convert to surface 
brightness units we adopt the FWHM beam sizes from | Aniano | 
|et al.| ( |2011| ) listed in Tab.[2] The (sub)mm data have also been 
converted to Jy/ D// assuming the appropriate effective beam sizes 
for the respective observations (see references in Table[3]). Since 
the standard calibration is done on point sources, but we are us- 
ing surface brightnesses to compile spatially resolved SEDs, we 
have applied the recommended extended emission calibration 
corrections to the SPIRE data (see SPIRE Observer's Manual: 
HERSCHEL-DOC-0798, version 2.4, June 7, 2011). 
PACS and SPIRE map background subtraction: The 
Scanamorphos "galactic" option preserves the large spatial 
frequency modes of the bolometer signal in the Herschel 
scan-maps. However, the absolute flux level of the background 
and spatial structures more extended than the scan-maps 
cannot be recovered since the exact level of the strong thermal 
background from the only passively cooled mirrors Ml and M2 
of the observatory is unknown. Therefore, it is not clear how 
the extended, large-scale emission levels in the final Herschel 
maps are related to the absolute flux level of the background, 
which is composed of emission features larger than the map 
size, the general galactic background ISRF, the cosmic mi- 
crowave background (CMB), and other possible contributions, 
all varying differently with wavelength. Since uncertainties in 
additive flux contributions affect flux ratio measurements and 
the respective derivation of physical parameters from the SEDs, 
in particular at low flux levels, we subtract the background 
levels from all Herschel maps, thus making them consistent 
with the chopped (sub)mm data (although chopping is a much 
more aggressive filtering than that present in the Herschel 
maps). The consequences and uncertainties of this bac kgr ound 
removal on the data analysis are discussed in Sects.|4.2|and|6.1 



7 



R. Launhardt et al.: The structure of low-mass cloud cores 



We emphasize that this approach was feasible only because our 
sources were already initially selected to be relatively isolated 
on the sky (see Sect.[2]). 

To accomplish this re-zeroing, we adop t a method simila r 
to that applied to Spitzer MIPS images in |Stutz etaL] ( |2QQ9a| ). 
For each object we identify a 4' x 4' region in the PACS and 
SPIRE images that is relatively free from spatially varying emis- 
sion and appears "dark" relative to the globule extended emis- 
sion levels. We impose the additional requirement that this re- 
gion is in or near a region in our complementary molecular line 
maps which is relatively free of 12 CO(2-l) emission. We always 
use the same region in all five Herschel scan-maps. For each 
band, we then calculate the representative flux level, /dc, in the 
4' x 4' region by implementing an iterative Gaussian fitting and 
cr-clipping scheme to the pixel value distribution at each wave- 
length. We do not consider pixels below 2cr from the mean and 
fit the main peak of the pixel value distribution. The /dc value is 
then determined by iteratively fitting a Gaussian function to the 
histogram of pixel values, where at each iteration we include one 
more adjacent higher flux bin. The final adopted value of /dc is 
defined as the mean value of the best-fit Gaussian that has the 
minimum cr value for all iterations. In this way we exclude pix- 
els with higher flux which effectively broaden the distribution 
and would cause a bias in the /dc calculation. The / D c values 
are then subtracted from the corresponding Herschel maps at the 
respective wavelength (see Table|4]for the /dc and cr values and 
coordinates of the reference region). 

Image convolution: Finally, all maps are convolved to the beam 
of the SPIRE 500 yum map, using the azimuthally averaged 
Herschel convolution kernels provided by Aniano e t al.| ( [20TT] ). 
For the convolution of the (sub)mm maps, we use a Gaussian 
convolution kernel of width equal to the difference in quadrature 
between the effective FWHM of the (sub)mm observations and 
the SPIRE 500 jam FWHM of 36! 4. 

After re-gridding to a common Nyquist- sampled pixel grid, 
we thus compile pixel-by-pixel SEDs for the wavelength range 
from lOOyimi through 1.2 mm. Photometric color corrections for 
each pixel and each band are iteratively derived from the 100- 
160yum and 250-500 yum spectral indices and polynomial fits to 
color correction factors provided for PACS in Table 2 of Miiller 
|eTafl { |20Tl) and in Table 5.3 of the SPIRE Observers Manual, 
and are applied in the subsequent SED modeling only. 

The nominal point- source calibration uncertainty of the 
PACS photometers listed in the latest calibration notes is 3% 
at lOOyurn and 5% at 160yum. However, additional uncertainites 
that are hard to quantify are introduced by, e.g., the conversion to 
surface brightness units for extended emission, the beam convo- 
lution (imperfect kernels), and the color corrections. Likewise, 
the final recommended calibration uncertainty of 7% for SPIRE 
does not account for some of the uncertainties introduced in 
the post-processing. Recent systematic comparisons between the 
Herschel and Spitzer surface brightness calibrations showed that 
they agree to within about 10%, with intrinsic uncertainties on 
both side^] Therefore, we adopt for all Herschel data the conser- 
vative value of 15% for the relative calibration uncertainty in the 
final maps. This number is used for calculating the weights of 
the individual data points and with respect to the ground-based 
data in the subsequent fitting procedure (Sect.|4.2|). 



4.2. Deriving dust temperature and column density maps 

The subtraction (or chopping out) of a flat background level from 
the emission maps implies that the remaining emission in the 
map at each image pixel is given by 

S v (v) = £1 (l - e" T(v) ) (B v (v, r d ) - 7 bg (v)) (1) 
with 

M d 

t(v) = N n m n —- K d (v) (2) 

M H 

where S v (v) is the observed flux density at frequency v, Q 
the solid angle from which the flux arises (here normalized to 
1 arcsec 2 ), r(v) the optical depth through the cloud, B v (v, T&) the 
Planck function, 7^ the dust temperature, Zbg(v) the background 
flux level, iV H = 2x N(U 2 ) + N(U) the total hydrogen column 
density, mn the proton mass, M d /M H the dust-to-hydrogen mass 
ratio, and K&(y) the dust mass absorption coefficient. and k& 
may vary along the LoS, in which case Eq.[T]has to be written 
in differential form and integrated along the LoS. However, for 
this overview paper we assume both parameters to be constant 
along the LoS and d iscus s the uncertainties and limitations of 
this approach in Sect. |6.2| 

For Afd(v), we assume for all sources in this paper the 
tabulated values listed by [Ossenkopf & Henning (1994]) for 
mildly coagulated (10 5 yrs coagulation time at gas density 
10 6 cm -3 ) composite dust grains with thin ice mantles (column 5 
in their Table 1, usually called "OH5"), logarithmically inter- 
polated to the respective wavelength where necessary. For the 
hydrogen-to-dust mass ratio in the Solar ne ighborhood, we adopt 
M H /M d = 110 (e.g., |Sodroski et"aT1|1997|). Note that the total 



gas-to-dust mass ratio, accounting for helium and heavy ele- 
ments, is about 1.36 times higher, i.e., M g /M^ « 150. 

The exact background flux levels at the location of the in- 
dividual clouds, /bg(v), are a priori unknown as explained in 
Sect. 4.1 At the frequencies relevant for this paper, this back- 
ground radiation is composed mainly of the CMB, the ex- 
tragalactic cosmic infrared background (CIB), and the diffuse 
galactic background (DGB). While the first contribution (CMB) 
is well-known and dominates at wavelength > 500 /mi, the mean 
levels of the CIB have been compiled by, e.g 



Hauser & Dwek 



7 see PICC-NHSC-TR-034 and PICC-NHSC-TN-029 under 
https://nhscsci.ipac.caltech.edu/. 



(2001 ) from flux measurements in the "Lockman Hole" by vari- 
ous space observatories (e.g., COBE and ISO) and ground-based 
(sub)mm instruments (e.g., SCUBA). However, at wavelengths 
< 350 yum, the DGB, which strongly varies with position, d omi- 
nates over C MB and CIB. Here we follow the ap proach of |Stutz| 
eTaL] pOlO) and use the |Schlegel eTaL] ( [T998] ) 100 yt/m IRAS 
maps and the ISO Serendipity Survey observations at 170yum to 
extrapolate approximate mean flux levels at the typical galactic 
locations of our sources. Since the uncertainty in the resulting 
temperature and column density estimates introduced by the un- 
certainty in the exact knowledg e of t he 4g levels is negligible, as 
discussed quantitatively in Sect. |6.1| we do not attempt to derive 
the local values of the DGB for each source separately. Instead, 
we adopt the following mean values for all sources: 0.3, 0.8, 0.5, 
0.3, 0.2, 2.9, and 6.5 mJy/arcsec 2 at A 100, 160, 250, 350, 500, 
850, and HOOyum, respectively. 

From the calibrated and background- subtracted dust emis- 
sion maps, we extract for each image pixel the SED with up 
to 8 data points between lOO/mi and 1.2 mm. These individ- 
ual SEDs are independently fit (x 2 minimization) with a single- 
temperature modified blackbody of the form of Eqs.[T]and|2]with 
r d and A^h being the free parameters. The individual flux data 
points are weighted with cr -2 , where cr is the quadratic sum of 
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flux times the relative calibration uncertainty (see Sect. |4.1[ ) and 
the mean rms noise in the respective map measured in regions 
of zero or lowest emission outside the sources (see Table]?]). In 
order to minimize the effect of T - Nu degeneracies in the fit- 
ting, because at the low temperatures in these cores, T < 10 K, 
even the 1 .2 mm emission does not represent the Rayleigh- Jeans 
regime where the SED slope is independent of T, the weight 
of data points with S v < cr was set to zero, i.e., these points 
were not considered in the fitting. To further avoid erroneous ex- 
trapolations of the SED fit to unconstrained shorter wavelengths, 
only pixels with valid fluxes in all 5 Herschel bands were fitted. 
This latter criterion was usually constrained by the PACS 100 yum 
maps. The best-fitting parameters were derived by employing a 
least squares fit to all flux values between lOOyum and 1.2 mm at 
a given image pixel, using a "robust" combination of the classi- 
cal simplex amoeba search and a modified Levenberg Marquardt 
method with adaptive steps, as implemented in the 'mfif tool of 
GILDAS0 

This procedure yields LoS-averaged dust temperature and 
colu mn d ensity maps of the sources, which are presented in 
Sect. 52_ The temperature maps provide a robust estimate of 
the actual dust temperature of the envelope in the projected 
outer regions, where the emission is optically thin at all wave- 
lengths and LoS temperature gradients are negligible. Toward 
the source centers, where cooling and shielding or embedded 
heating sources can produce significant LoS temperature gradi- 
ents and the observed SEDs are therefore broader than single- 
temperature SEDs, the central dust temperatures will be over- 
estimated (in the case of a positive gradient in cold sources) or 
underestimated (negative gradient in internally heated sources). 
For these reasons, the column density maps, which come out 
usually very smooth, are corrupted and exhibit artifacts in re- 
gions of ~20" -40" radius around the warm protostars (see, e.g.. 
Figs. 30] or 
Sect. 



29). These caveats are discussed and quantified in 



6.2 including some already worked-on solutions. 



5. Results 

5. 1. Herschel FIR maps 

The resulting calibrated dust emission maps at X 100, 160, 250, 
350, and 500 yum and at original angular resolution are presented 
in Figs. [12] through [23] The Herschel maps are accompanied 



by optical (red) images from the second Digitized Sky Survey 
(DSS20), which were obtained through the Sky Vie\^\ interface . 
These optical images clearly show the regions of highest extinc- 
tion as well as, in several cases, extended cloudshine structures. 
All maps are overlaid with contours of the (sub)mm dust contin- 
uum emission observed with ground-based telescopes. 

The 160yum through 500 yum maps are usually very similar 
to each other in appearance, outlining the thermal dust emission 
from the dense cores along with the often filamentary or tail- 
like, more tenuous envelopes. Except for some filamentary or 



8 http://www.iram.fr/IRAMFR/GILDAS 

9 The Digitized Sky Surveys were produced at the Space Telescope 
Science Institute under U.S. Government grant NAG W-2166. The im- 
ages of these surveys are based on photographic data obtained using the 
Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt 
Telescope. The plates were processed into the present compressed dig- 
ital form with the permission of these institutions. 

10 SkyView has been developed with generous support from the NASA 
AISR and ADP programs (P.I. Thomas A. McGlynn) under the auspices 
of the High Energy Astrophysics Science Archive Research Center 
(HEASARC) at the NASA/ GSFC Astrophysics Science Division. 



tail-like extensions like, e.g., in CB 17 (Fig.[l4]), the -18' SPIRE 
maps usually cover the entire extent of the detectable FIR dust 
emission down to the more or less flat background levels. The 
slightly smaller PACS maps (-10') also usually cover the entire 
extent of the clouds, but in some cases cut off a bit more of some 
filamentary or tail-like extensions (e.g., Fig.[T4|). 

The PACS lOOyurn maps often show significantly fainter or 
less extended emission than the maps at longer wavelengths, ow- 
ing to the fact that lOO/mi samples the steep short- wavelength 
side of the SED at the low temperature of the dust in these clouds 
(Fig.[4]). In some cases, the diffuse 100 yum emission even ex- 
hibits a "hole" or "shadow" at the location of the column density 
peak (e.g., CB 27; Fig. 16), hinting at extremely cold and dense 
(high column density) cores with thin warmer envelopes, which 
are thus go od candidates of cores at the verge of protostellar col- 
lapse (see |Stutz et al.|2009a| for a discussion of 24 yum and 70 yum 
shadows in dense cores). 

On the other hand, the lOOyurn maps, which also have the 
highest angular resolution of all our Herschel maps (« 7"), are 
most sensitive to embedded heating sources like, e.g., protostars. 
Although even the previously known very low luminosity object 
(VeLLO) in CB 130 ( |Kim et al.|[20TT| see also |Dunham et al.| 
2008 ) is well-detected (Fig.[20]), our maps do not reveal any ob- 
vious previously unknown warm compact source in any of the 
globules, confirming that all our presumably starless cores are 
indeed starless. The only possible exception might be CB 17 (see 
Fig.[T4|), where we find hints of an extremely low-luminosity 
(Lboi < 0.04 L Q ) cold embedded source very close (in projec- 
tion) to CB17-IRS ( |Chen et al.|20l2||Schmalzl et al.|2013| ). 

5.2. Dust temperature and column density maps and 
radial profiles 

The LoS-averaged dust temperature and hydrogen column 
density maps derived with the fitting procedure described in 
Sect. |4.2| are presented in Figs. [24] through [35] All sources 
show systematic dust temperature gradients with cold interi- 
ors (1 1 - 14 K) and significantly warmer outer rims (14 - 20 K). 
Embedded heating sources like protostars, including the VeLLO 
in CB 130 (Fig. [55} Kim et al.|2011| ), show up very clearly in the 
temperature maps. 

To characterize the column density profiles quantitatively, 
we fit circular profiles around the col umn densi ty peak to the 
maps, using a "Plummer-like" profile ([Plummer 191 1 1 see also 
Whitwo rth & Ward-Th ompson 2001}, modified by a constant 
term to account for the obse rved outer column density "plateau" 
(see discussion in Sect.|5.7|): 



AN 



Nn (r) = (mr/nfY 12 



+ N mtt if r < r n 



if r > r out 



The peak column density is then 
N = N n (r = 0) = AN + N out . 



(3) 



(4) 



This profile accounts for an inner flat (column) density core at 
r < r\ where 



AN 

Ni =Nu(n) = ^ +AU 



(5) 



approaches a power-law with index p at r » r\, turns over into 
a flat outer column density "plateau" outside 



l(AN\ 2/p 



(6) 
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Table 5. Column density and dust temperature profiles 



Source 


N b 


iV out 


n b 


r 2 b 


r b 
• out 


P b 


r ~p b 
1 in 


r T 1 C 

1 peak 


^ out 




r 9t 

[cm 2 ] 


r 9t 

[cm 2 ] 


[AU] 


[AU] 


[AU] 




[K] 


[K] 


[K] 


LB 4 


7.5E21 


0.UL19 


1 CC A 


o.Zh4 


C\ CC A 

9.jL4 


J.U 


1 A A 

14. U 




1 A C 

19. j 


o 1 n OA TAT 

LB 1 / - oMM 


z.jbzz 


1 nco 1 

l.Ubzl 


1 AC /I 


1 1 11? A 

3.1L4 


Q C C A 


J.U 


1 1 T 
11. / 




1 Q 1 

1 J.l 


CB zo - bJVlJVlz 


1 .4ilzz 


1 OC0 1 


5.4iiJ 


O 1C /I 

z.zh,4 


3.4ii4 


O A 

ZA 


1 o o 

Iz.o 




1 /I A 

14. U 


Cr$ z/ 


i.yc,ZZ 


z. ibzi 


1 A.UA 

l.or!#4 


1 1 C/l 

z. 1H4 


1 Ot? /I 

Z.OH4 


a n 
4.U 


1 1 o 

11. y 




1 A 1 

14. z 


O /CO 

r$ Do 




4.UiiZU 


1.UH4 


Z. lrL4 


/I 


j.U 


1 1 Q 

YV.y 




1 o c 
1 O.J 


Z44 - olVllvlZ 


4.0r!/ZZ 


1 . 1HZ1 




O.Ui}4 


Q HTh/I 

>>.Ui14 


l.o 


1 1. J 




1 /I ^ 
14. J 


CB130-SMM J 


2.6E22 


2.0E21 


7.2E3 


2.8E4 


4.7E4 


1.8 


12.1 




13.5 


Cr> ' J 


j.Uiizl 


j.UiizU 






1.4r!/j 




1 A ^ 
14. j 


1 Q /I 

lo.4 


1 J. / 


ZD - olVllvl 1 


1 fYE"")'") 

1 .UUZZ 


q ninOn 
y.UCfZU 












ion 


1 /I £ 
14.0 


BHR 12^ 


5.4E22 


1.0E21 


1.4E4 


J). 


1.1E5 


4 


14 8 

IH-.O 


1 O.J 


1 7 ^ 


CB 68 


2.4E22 


5.0E20 


1.1E4 


2.7E4 


2.9E4 


4.0 


12.7 


19.0 


16.5 


B335 


3.1E22 


6.2E20 


3.6E3 


1.8E4 


3.8E4 


2.3 


14.7 


18.8 


16.2 


CB 230-SM1VK 


2.7E22 


5.0E20 






1.4E5 




14.5 


19.6 


15.6 


CB 244-SMM1 


1.5E22 


1.0E21 


1.2E4 


5.5E4 


9.0E4 


1.7 


11.5 


19.6 


14.2 


starless cores 


2.3+1.0E22 


1.1+0.6E21 


1.2+0.6E4 


3.5+1.6E4 


5+2E4 


3.3+1.2 


12.2+0.8 




15.3+2.2 


protostellar cores 


2.4+1.4E22 


0.7+0.2E21 


0.9+0.5E4 


3.4+1.2E4 


9+4E4 


3.0+0.9 


13.6+1.2 


19.0+0.4 


15.7+0.9 


all 


2.3+1.3E22 


0.9+0.5E21 


1.1+0.5E4 


3.5+1.5E4 


7+4E4 


3.2+1.2 


13.0+1.2 




15.5+1.8 



Notes. In contrast to the other tables, sources in this table have been grouped into st arles s and protostellar. {a) Coordinates are listed in Table[6] 
ib) Derived from circular fits to the dust temperature and column density maps; see Sect. 52_ Eqs.|3]through[9] (c) Peak temperature at the position of 
the embedded heatig source (protostar), severely underestimates the true central dust temperature because of beam smoothing and LoS averaging. 
(J) CB 130: the embeded VeLLO causes only a very small temperature increase in the smoothed dust temperature map. (e) Bad fit due to large 
unresolved LoS gradients, no meaningful values for some of the parameters. (f) The embedded double sources in BHR 12 and CB 230 are not 
resolved in the smoothed dust temperature maps. 



where 

N 2 = N u (r 2 ) = 2xN 0] 



(7) 



and is cut off at r out . This outer boundary of the thin envelope or 
halo is not well-recovered in most cases because filamentary ex- 
tensions emphasize deviations from circular / spherical geometry 
with increasing distance from the core center and the flux levels 
gradually decline below the noise cut-off (see Sect. |4.2| ). For the 
same reason, the flatness of the outer column density profile at 
level A^out might be partially an artifact of circularly averaging 
noncircular structure at a level just above the noise. 

To empirically fit and parameterize the radial temperature 
profiles of purely externally heated cores, we use a similar profile 
of the form: 



T(r) = T out - 



AT 



v2\*/ 2 ' 



(8) 



The central temperature minimum is then given by 

T in = T 0Ut -AT . (9) 

Before fitting radial profiles to the maps, we mask some of 
the spurious low-SNR edge features and tail-like (asymmetric) 
extensions, always applying the same spatial mask in both the 
column density and dust temperature maps. For cores with em- 
bedded heating sources (protostars), we also mask the region 
with local temperature increase before fitting the radial profiles 
with Eqs.[3]and[8]and also list in Table[5]the value of the central 
temperature maximum r peak . The radius of the masked region 
varied between 20" and 40". This extrapolation for A^o and T m 
of the outer profile into the core center, where both the local tem- 
perature and the column density estimates are very uncertain due 
to large and unresolved LoS temperature gradients, leads to an 
increased uncertainty in particular of the value for in the pro- 
tostellar cores. This and other shortcomings of this model and 
the interpretation of its parameters are discussed in Sect.[rjT2 



To illustrate the quality of these radial profile fits and the 
meaning of the various parameters, we show in Fig. [3] the radial 
distribution of A/h and T d values in the resulting maps for B 68 
along with the best-fit profiles and the parameters labeled on the 
diagram axes. The best-fit parameter values for A^, A^ out , r\, r 2 , 
r out , p, r^, r peak (for cores with embedded heating sources), and 
r ou t of all sources are listed in Table[5] We do not list parameters 
r T and q, but note that from the core centers outward, the LoS- 
averaged dust temperature typically raises by 1 K at a radius of 
(1 +0.3) 10 4 AU. 

Note that in the next section, we define core and cloud sizes 
via the mean extent of certain column density contours in the 
maps, and not via the radial profile fit parameters r\ and r out . In 
particular for very elliptical sources and sources with extended 
halos, the derived mean diameters are therefore not necessarily 
identical to twice the corresponding radial profile-fit radii. 



5.3. Source morphologies, integral properties, and SEDs 

In addition to the pixel-by-pixel SEDs and the column density 
and dust temperature maps, we also derive integrated source pro- 
perties using certain column density levels to define the integra- 
tion areas. Total cloud sizes are derived by measuring the major 
and minor extent of the Af out contour in the column density maps, 
which is marked as the outer yellow contour in Figs. [24] through 
35 The resulting linear diameters (invoking the distances listed 
in Tab.[T]) are listed in Tab. [6] The mean radii of the globules in 
our sample range from 0.1 pc (CB 130) to 0.5 pc (CB 230), with 
a mean of 0.22 + 0.1 pc (4.8 + 2.2 10 4 AU), which is about twice 
the typical Jeans lengths of these clouds (~ 0.13 pc, assuming 
the mean radius and mass and a me an temperature of 13 K; see 
also stability discussion in Sect. ^5 ). The mean projected aspect 
ratio of the clouds is 0.66 + 0.16, with the most extreme cases 
CB 6 and BHR 12 (0.4) and B 335 (0.98). These measures do 
not account for the long tails, e.g., in CB 17 or CB 68, nor do 
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Table 6. Physical parameters of globules and embedded cores 



Source 


R.A.,Dec. (J2000) fl 


Size^ 


b:a^ 


M H C 






^smm C 
Lbo\ 


Morphology^ 


Fig 


ures 








J 


[AU] 




[M Q ] 


[L Q ] 


[K] 


[%] 


Evol. class 








CB4 






1.2E5 


0.74 


1.6 


2.8 






R, T(0.5) 


12 


A.l 


24 


CB 4 - SMM 


00:39:05.2, 


+52:51:47 


5.2E4 


0.62 


0.7 


0.77 


23 


11.7 


starless, unbound 








CB6 






1.2E5 


0.40 


5.9 


3.5 






C, T(2.6) 


13 


A.2 


25 




CB 6 - SMM 


00:49:24.3, 


+50:44:36 


3.9E4 


1.0 


1.0 


2.5 


158 


5.1 


CI. I 










CB17 






8.6E4 


0.67 


3.0 


1.3 






C, T(>1.0) 


14 


A.3 


26 




CB 17- SMM 


04:04:37.7, 


+56:55:59 


2.1E4 


0.86 


1.0 


0.37 


18 


22.3 


starless, stable 










CB17-IRS 


04:04:34.0, 


+56:56:16 








0.12 


80 


2.2 


CLP 








CB 26 






7.3E4 


0.75 


2.4 


1.9 






C, M, T(0.6) 


15 


A.4 


27 




CB 26- SMM 1 


04:59:49.3, 


+52:04:39 


2.0E4 / 


1.0 / 


0.3' 


0.47' 


81' 


6.5' 


CI. I 










CB 26-SMM2 


05:00:14.5, 


+52:05:59 


2.2E4 / 


1.0' 


0.6' 


0.30' 


20' 


17.3' 


starless, stable 








CB 27 






6.7E4 


0.77 


3.7 


2.2 






C, T(0.3) 


16 A.5 


28 




CB 27 -SMM 


05:04:08.1, 


+32:43:30 


2.6E4 


0.61 


1.2 


0.44 


18 


20.5 


starless, stable 










BHR12 






1.2E5 


0.39 


11 


29 






C, M, T(0.9) 


17 


A.6 


29 




BHR12-SMM 


08:09:32.7, 


-36:05:19 


3.1E4S 


0.82* 


3.7* 


14.9^ 


104* 


3.2^ 


C1.0/I 7 * 










CB 68 






6.5E4 


0.55 


2.1 


2.1 






C, T(0.7) 


18 


A.7 


30 




CB 68 - SMM 


16:57:19.7, 


-16:09:23 


1.5E4 


0.79 


0.5 


0.86 


41 


6.3 


C1.0 










B68 






5.6E4 


0.83 


1.3 


1.3 






c 


19 


A.8 


31 




B 68 - SMM 


17:22:38.3, 


-23:49:51 


1.6E4 


0.88 


0.6 


0.23 


19 


17.4 


starless, stable 










CB130 






4.9E4 


0.52 


2.7 


1.6 






M 


20 


A.9 


32 




CB 130 -SMM 


18:16:15.6, 


-02:32:45 


2.0E4 


0.63 


1.0 


0.37 


38 


20.8 


VeLLO* 










B335 






4.8E4 


0.98 


0.8 


1.3 




4.4 


R 


21 


A.1C 


> 33 


B 335 - SMM 


19:37:00.7, 


+07:34:08 


8.5E3 


0.80 


0.2 


0.57 


37 


C1.0 








CB 230 






2.1E5 


0.69 


18 


25 






C 


22 


A.ll 


34 


CB 230 -SMM 


21:17:39.9, 


+68:17:36 


3.5E4 


0.55 


3.0 


10.4 


189 


3.7 


CI. I* 








CB 244 






1.4E5 


0.61 


12 


8.2 






M 


23 


A.12 


35 


CB 244- SMM 1 


23:25:47.3, 


+74:17:44 


2.0E4 Z 


0.87 z 


0.8 Z 


1.82' 


56* 


4.3 Z 


C1.0 








CB 244-SMM2 


23:25:26.8, 


+74:18:22 


2.3E4 


0.87 


2.1 


0.53 


17 


24.2 


prestellar 









Notes. (a) Coordinates of column density peak (starless cores) or dust temperature peak (embedded source). (b) Mean diameter ( V« x b) and aspect 
ratio of the Af out column density contour listed in Table[5]for the globules (i.e., « 2 X r out , Eq.[3]and Fig.[3]), and of the Ny contour (Eq.[l0) for the 
embedded cores (i.e., 2 X ry, Eq.[TTJ. (c) Measured inside the Af out contour (globules) or the Ny contour (cores), except where marked otherwise. 
(J) Morphology description based on visual inspection of DSS2-red optical images and Herschel dust emission maps: R: round or slightly elliptical; 
C: cometary- shaped; M: multiple cores; T(x): tail of length xpc. (e) CB 17 - IRS is identified in the PACS 100 //m map and as a local temperature 
maximum in the T-map, but cannot be separated from SMM i n the column density map. ( ' } CB 26 - SMM1 and SMM2: ellipses, marked by dotted 



lines in Fig. 27 were use d ins tead of the Ny contour (cf. Sect. |5.3| >. (g) BHR 12: An ellipse, marked by dotted lines in Fig.|2 9j was used instead of 
the Ny contour (cf. Sect.[53]). {h) BHR 12: The two embedded sources SMM1 (Class I) and SMM2 (Class 0; sep. 21"; see Laun hardt et al.|2010| ) 
are not resolved here. (z) CB 130: additional low- mass prestellar core ~3 0" west and Class I YSO ~15" east of Class core. w CB 230: the two 
embedded sources IRS1 and IRS2 (sep. 10"; see Launhardt et al. |2010} are not resolved here. (/) CB 244 - SMM1: an ellipse, marked by dotted 
lines in Fig. [35] was used instead of the Ny contour (cf. Sect. |5.3| ». 



they accurately reflect the actual aspect ratios of the individual 
sources, since the 3-D projection of the elliptical clouds onto the 
plane of sky remains unknown. Note that these and the following 
uncertainty values only represent the variance of the mean and 
do not include systematic uncertainties, which are discussed in 
Sect.0 

Since cometary tails in gl obules appear to be failry ubiqui- 
tous (see Tab.[6]and Sect. |5.8| ), we have also assessed the possi- 
bility that some of the apparently tail-less globules (e.g., B 335 
or CB 130) could possibly also have a tail pointing at us or away 
from us. However, given the length of the tails observed, e.g., 
in CB 6 (Fig.|A^]), CB 17 (Fig.[X3]), or CB 68 (Fig.[A77]), it is 
extremely unlikely that a globule like B 335 has a tail that is di- 
rectly aligned with the line of sight. Another argument against 
this possibility comes from the fact that, at least in the starless 
cores, the (LoS-averaged) temperature minimum always agrees 
with the column density maximum. Should a long tail, like that 
seen in e.g., CB 17, be oriented along the line of sight, then its 
total column density would sum up to a value similar to or even 



larger than that observed toward the core centers. However, since 
the tails are not as well-shielded from the IRSF as spherical cores 
and therefore have approximately the same average temperature 
as the envelopes, such a projection could not result in a clear 
temperature minimum at the position of the column density max- 
imum, as observed in all starless cores. 

Total hydrogen masses are derived by integrating the col- 
umn density maps within the N out contour and are also listed 
in Tab.[6] These masses range from 0.8 M Q (B 335) to 18M 
(CB 230), with a mean value of <M c i oud ) = 5.4 + 5.0 M . Using 
the somewhat more conservative column density threshold N2 
(Eq.[7]) would have the benefit of being less susceptible to inho- 
mogeneities of the background level and having closed contours 
in all sources, but would cut off some of the extended emission 
that most likely originates from the globule. Using N2 instead of 
N out would lead to -10 - 15% smaller masses and sizes for the 
globules compared to the values listed in Tab.[6] 

Integral cloud SEDs (not shown) were obtained by integrat- 
ing the various emission maps within the N out contour. Where 
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1300 



A [/-im] 

100 12 



2.2 



1300 



A [/xm] 

100 12 



2.2 



1300 



A [/U.m] 
100 12 2.2 



-12 



-16 



-12 



-16 



-12 



-12 



■— 1 -14 



-12 



-16 



CB4 

/ 1 

" / 1 

_ / T bol - 23 K 

star-less 
i i i 


CB17-SMM 

; A ; 

_ ' T bol = 18 K 

star— less 
i i i 


CB26-SMM2 

T bol = 20K 
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Fig. 4. Integrated spectral energy distributions of the dense cores and embedded sources listed in Table^ Fluxes are integrated 
within the N\> contour (Eq. 10), as are the core masses and luminosities listed in Table[6] Solid lines between 3 mm and lOOyurn 
mark the modified blackbody fits to the data points, while dashed lines only represent a logarithmic interpolation of the data points 
shortward of lOO/mi (see Sect.|5.3|). 



available and appropriate, we also included ground-based op- 
tical and NIR data, Spitz er IRAC and MIPS m aps, as well as 
ISO and IRAS fluxes (see |Launhardt et al.|2QT0] ). The resulting 
SEDs where fitted between A 100//m and 1.2 mm in the same 
way as the individual image pixel SEDs (Sect. |4.2| ). To facil- 
itate integration of the SEDs at wavelengths shortward of the 
SED peak, fluxes at X < 100//m were logar ithmically inter- 
polated as described in Launhard t et al.| ( [2010| ). Bolometric lu- 
minosities (Lboi) were derived by integrating the SEDs over the 
complete wavelength range. Submillimeter luminosities (L smm ) 



were derived by integrating the SEDs at wavelengths longward 
of 350ytmi. The resulting total cloud luminosities range from 
1.3 L Q (CB 17, B 68, B 335) to 29 L Q (BHR 12), with a mean of 
(^cioud) ~ 6.7L . See Tableland Fig. [7] for the distribution of 
masses and luminosities. 



In four out of the twelve globules, the Herschel and submm 
maps reveal multiple cores (CB 26, BHR 12, CB 130, CB 244 - 
all known before to be multiple), while the other eight globules 
appear single-cored at the Herschel resolution. In the tempera- 
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Fig. 3. Radial profiles of column density and LoS-averaged dust 
temperature of B 68 (cf. Fig. 9 of |Nielbock lFalJ|20T2] for the 
ray-tracing results of the same data). Small black dots show the 
data pixel values over radial distance from the column density 
peak (see Fig. [3T|) Solid gray lines show best fits to these data 
with Eqs.[3]and[8jand the parameters listed in Table[5] The pa- 
rameters are also labeled on the diagram axes to illustrate their 
meaning, with the following characteristic radii: r\ = flat core 
profile radius (Eqj3l, ry = radius to define consistent dense core 



boundary (Eqs.|10|and 1 1 ), r 2 = cloud radius at transition to halo 
(Eq.[6]), and r out = outer cloud boundary (Eq.[3]). 



ture and column density maps, the double cores are only resolved 
clearly in CB 26 and CB 244. 

The identification of a common criterion for characterizing 
the extent and measuring masses and luminosities of the dense 
cores was less straight-forward. Since the (column) density pro- 
file of a single-cored globule continuously decreases from the 
center toward the outer boundary of the globule, there is no 
simple way to observationally distinguish between a "core" and 
an "envelope", like in larger molecular clouds where the outer 
boundary of a core is usually judged to be where the (column) 
density profile merges with the surrounding cloud level. The 
flat core radius, r\, of the column density profile is also not a 
good criterion since its value depends on the fitted power-law in- 
dex p (Eq.[3]), which is not well-constrained in most cases since 
V2 ^> ri, i.e., the range where the column density profile ap- 
proaches a power-law and p can be derived is too small in many 
cases. Therefore, we define an empirical, but well-reproducible 
column density threshold for all dense cores by 



AN 

Ny =Nn(n>) = — +Af 0Ut , 
e 

where e is Euler's number and 
r v = n ^e 2 'P - 1 . 



(10) 



(11) 



(Eq.[5|. This column density threshold is high enough to per- 
mit a reasonable separation of the subcores in CB 26 (Fig. [27]) 
and CB 244 (Fig. [35]), while at the same time being low enough 
to fully include the regions of corrupted column density values 
around embedded protostars (due to strong unresolved LoS tem- 
perature gradients; cf. Sect. |4.2| ). Note that this is an observa- 
tionally driven definition of "core" that is not based on a well- 
characterized physical transition between core and envelope (see 
also discussion in Sect. |5.4| ). In a few cases, we had to use el- 
lipses to define the core areas, which were however guided by 
the Ny contour. In CB 26 (Fig. 27 ), the Ny contour still includes 
a bridge between cores SMM1 and SMM2, such that an artifical 
cut-off by an ellipse was necessary. In addition, the corruption of 
the column density map around SMM1 by unresolved LoS tem- 
perature gradients leads to an offset between the column density 
peak (artifact) and the temperature peak (actual location of the 
embedded YSO). This latter effect required to also use ellipses 
for CB 244-SMM1 (Fig.|35) and BHR 12 (Fig.|29j. 

The resulting integral SEDs of the core regions enclosed 
by the Ny contours or the respective ellipses (where marked in 
Figs. 24 through 35 ) are shown in Fig. [4] These SEDs where fit 
in the same way as described above for the integral cloud SEDs. 
The integral properties derived from these subregions are much 
more representative of the cold starless cores or warm proto- 
stellar cores than the total globule quantities which include the 
extended envelopes that are dominated by external heating from 
the ISRF. Resulting sizes, aspect ratios, masses, and bolomet- 
ric luminosities of the cores are listed in Table[6] The relation 



For p = 2/ln(2) (-2.9, which is actually very close to the 
mean (p) of our sampe), we obtain Ny = N\ and ry = r\ 



between SEDs and evolutionary stages is discussed in Sect.[54 
The mean core diameters, as defined by the Ny contour and 
listed in Table^ range from 8.5 x 10 3 AU (0.04 pc; B 335) to 
5.2 x 10 4 AU (0.25 pc; CB 4), with a mean of (2.5 ± 1) x 10 4 AU 
(0.12 ± 0.05 pc). The mean aspect ratio of the cores is (b : a) = 
0.8 ± 0.14, i.e., they are somewhat rounder than the globule en- 
velopes. However, this latter difference is not significant since 
beam- smoothing does affect the cores more than the envelopes. 
The mean angular diameter of the cores derived this way is for all 
sources more than twice larger than the beam (HPBW 36'.'4), and 
for half of the sources more than three times larger. Hence, the 
effect of the beam size on the core definition is small, although 
not completely negligible in some sources. A more severe lim- 
itation on the accuracy of core mass estimates may come from 
the fact that in the presence of unresolved warm protostars and 
large LoS temperature gradients, the column density maps are 
corrupted at the location of the protostar and the peak column 
density can be underestimated (see Sect. |4.2| ). In these cases, the 
core masses are underestimated and the core sizes have larger 
uncertainties. 

The mean total hydrogen mass of all dense cores in our sam- 
ple is 1.2 ± 0.9 M , with the lowest-mass core being in B 335 
(0.2 M©), the highest-mass core in BHR 12 (3.7 M©), and no 
significant systematic difference between starless and protostel- 
lar cores. The cores in single-core globules contain on average 
25+6 % of the total mass of the globule, irrespective of whether 
they already contain a protostar or not. In globules with two 
cores (CB 26 and CB 244), the total core mass fraction is compa- 
rable (30+5 %), but the individual cores constitute only 15+6 % 
of the total globule mass. In terms of their luminosity, starless 
and protostellar cores are more distinct from each other. The 
mean bolometric luminosity of the seven starless cores (Tab.[6] 
including CB 130) is 0.43 + 0.15 L , while that of the cores with 
embedded protostars is 4.5 L© with a large scatter ranging from 
0.5 L to 15 L (excluding CB 17 - IRS). The dense starless cores 
emit on average 20+5 % of the total luminosity of their host- 
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ing globules, while protostellar cores emit on average 40+10 %. 
These relations are discussed in more detail and in the context 
of the thermal structure of the globules and the relative effects of 



external and internal heating in Sect. 5.6 



5.4. Evolutionary stage tracers 

The spectral classes and the approximate evolutionary stages of 
the sources in our sample were already known from earlier ob- 
servations and modeling (see Tab.p]and |Launhardt et al.|2Q10j ). 
With the Herschel data in hand, which provide a much more 
complete coverage of the thermal SEDs and allow for more ro- 
bust mass, luminosity, and r bo i estimates than befor e, we re- 
ev aluate the classic al evolutionary tracers L smm / Lb i ( Andre et 
aL][T993]) and r bo i ( [Myers & Laddl[T993] ). Figure^ shows the 
L S mm / Lboi ratio vs. Tboi for all dense cores in the 12 globules of 
our sample. Uncertainties of these values are not straigt-forward 
to assess because the formal error bars are much smaller than 
the effect of exactly how a core and the respective integration 
area are defined. Based on tests with slightly different core def- 
initions we estimate the relative uncertainties to be about 10%. 
In contrast to earlier "pve-Herscher papers, the new L smm / L bo i 
vs. Tboi diagram now also includes purely externally heated star- 
less cores. At the end of this section, we also use complemen- 
tary information to relate the empirical "Class" of the individual 
sources to their actual evolutionary stage. 
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Fig. 5. Lboi/^submm ratio vs. bolometric temperature of the glob- 
ule cores (see Tableland Sect. 54). The size of the symbols 
scales with the peak column density of the respective sources. 
Error bars on CB 244 - SMM1 illustrate the 10% relative uncer- 
tainly on both values (Sect. |5.4| ). The "classical" Tboi boundaries 
for Class protostars are indicated by the vertical dashed lines. 
The two triangles, connected by a dotted line, represent one syn- 
thetic protostar - envelope system taken from (Robitaille et al. 
2006}, once seen close to edge-on (Tboi = 34 K) and once seen 
nearly pole-on into the outflow cavity (Tboi = 119 K; see discus- 
sion in Sect.|5.4b. 



All starless cores in our sample, i.e., cores that do not have a 
compact lOOyurn source or any other signs of central heating or 
star formation, have Tboi < 25 K and 10% < L smm / Lb i < 30%. 



We also find a correlation between Tboi and L smm / Lb i ratios 
and both mass and peak column density of starless cores in the 
sense that colder cores (lower r bol ) with higher L smm / L bo i ra- 
tios tend to be more massive (see Fig. [5}. Since these cores have 
no significant internal heating sources, this correlation is likely 
to reflect purely the degree of external heating by the ISRF and 
shieldin g ra ther than an evolutionary effect (see also discussion 
in Sect. 5.7). Massive cores with higher column densities are 
better shielded and can cool down more than less-massive and 
less- shielded cores. Even in the unlikely case that some of these 
cores would already have undetected first hydrostactic cores 
(FHSCs), their existence would have no measurable effect on 
Tboi or the L smm / Lb i ratio. Furthermore, an evolution of these 
integral quantities through the lifetmine of the FHSC is for the 
same reaso ns also not expected theoretically (e.g., Commercon 
|et al.|2012| ). On the other hand, cores in larger molecular clouds 
with active star formation may be expected to have both more 
shielding and a stronger local ISRF, which may lead to different 
temperature profiles and different location in the L smm / Lb i vs. 
Tboi diagram as compared to isolated cores. 

All cores with embedded sources previously classified as 
Class protostars are confirmed to have 30 K < r bo i < 70 K and 
and have 3% < L smm / Lb i < 7%. The only exception is the core 
of CB 130 with the embedded VeLLO, which has L smm / Lb i ~ 
21%. All cores with embedded sources previously classified 
as Class I YSOs are confirmed to have Tboi > 70 K and have 
2%<L smm /L bo i<7%. 

While the L smm / Lb i ratios we derive for the isolated Class 
sources agree well with the uppe r half of range of the values 
listed by, e.g., |Andre et al.| ( [2000| ) for Class sources in larger 
molecular clouds, we find significantly higher (up to a factor of 
10) L smm /L bo i ratios for the more evolved Class I YSOs than 



derived and propos ed as threshold by Andr e et aL] ([1993 ) and 

< 0.5%). In fact, we do not 



Andre et al. 



(2000 



(Class I) 



smm / ^bol 

rind a significant systematic difference in the L smm / Lb i ratios 
between Class p rotostars and Class I Y SOs, but see at most a 
slight trend. As in Launh ardt et al.| ( [20T0] ), we argue here that this 
systematic difference in L smm / L^ \ ratios between isolated and 
embedded Class I YSOs is unlikely to be an artifact of distance 
bias, angular resolution, or SED coverage, but rather reflects a 
combination of different methods of identifying the core bound- 
aries and possibly real environmental differences. Luminosities 
for these earlier embedded samples were derived by |Andre & 
Mon tmerle| ( |1994| ), [Moriarty-Scnieven et aL] ( |1994| ), and others 



based on IRAS point source fluxes and ground-based 800 and 
HOOyum or 1.3 mm maps only. The compact, warm embedded 
sources dominate the small submm maps completely and the re- 
maining low-level emission from a possible envelope was most 
likely chopped out to a large degree against the extended emis- 
sion from the surrounding cloud (inter-core) material. Hence, 
these older observations could not correctly account for the en- 
velope emission. In contrast, the column density profile of a 
globule core continuously decreases toward the "edge" of the 
small cloud, the isolated envelopes are heated more by the ISRF, 
and we do have spatially resolved emission maps at wavelengths 
between lOOyurn and 1.3 mm at hand. The net effect is that we 
recover more envelope emission than those earlier observations 
and the envelopes of the isolated sources are more luminous than 
those of the embedded sources, thus elevating the L smm / L bo i ra- 
tios. 

When trying to interpret the quantities Tboi and L smm / Lb i as 
evolutionary tracers, one has to keep in mind that, as soon as an 
outflow starts to open a cavity in the envelope, i.e., already early 
in the Class collapse phase, the observed SED at wavelengths 
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shortward of the thermal peak at A « 160 yum become s highly de- 
pendent on the viewing angle toward that cavity (e.g., | Whitney et| 
|al.|2003||Enoch et al.|2009[|Offner et al.|2QT2] ). Correspondingly, 
as long as one does not take into account the full source ge- 
ometry, but instead simply integrates the observed fluxes over 
47T, quantities derived from the SED also become projection- 
dependent. To illustrate this, we have taken one specific syn- 
thetic protostar - envelope system from Robitaille et al. ( |2006 



model ID: 3008984) with properties similar to our "Class 0' 
sources (stellar and envelope mass, luminosity, envelope radius, 
etc.), and analyzed its synthetic SED in the same way as our 
real sources. This synthetic source was once "observed" close to 
edge-on (i - 87°, w.r.t. the outflow axis), and once close to pole- 
on (i = 18°). Although the absolute L smm / Lb i ratios derived for 
this model are systematically somewhat lower than those of our 
real sources, possibly due to a combination of the simplicity of 
the specific model and the same physical reasons as discussed 
above, the r bo i values for the two projections (34 K and 119 K) 
and the L smm / Lb i trend cover the full range of "Class 0" and 
"Class I" sources (see Fig.[5]). Hence, without involving more 
specific information on the opening angle and orientation of 
the outflow cavities, e.g., from NIR images and high-resolution 
molecular line maps, together with radiative transfer modeling, 
a robust evolutionary classification of sources with r bo i > 30 K 
just from the SEDs is not possible. Protostars seen pole-on have 
SEDs that can be very similar to those of more evolved YSOs 
("Class I" stage) seen edge-on, i.e., for individual sources there 
is no clear one-to-one correspondence between spectral Class 
and evolutionary stage. 

For CB 26- SMM1, we know that it is seen nearly edge- 
on ( [Launhardt et al.|2009| ), i.e., its classification as YSO in the 
Class I stage seems to be robust. BHR 12 is harder to judge since 
the two sub-cores with nonaligned outlows are not resolved in 
our SED. However, the NIR images also seem to indicate that we 
are viewing them from a more edge-on orientation (e.g., Hodapp 
|& La dd 1995). Hence, its classification as YSOs in the Class I 
stage seems also robust. For CB 6 and CB 230, the NIR images 
suggest that we are at least partially looking into the outflow 
cones (Launhar dt et al.|20 10), i.e., they could well be protostars 
seen pole-on. 

5.5. Stability of starless cores 

In order to assess which of the starless cores are unstable against 
gravitational collapse und must therefore be prestellar in nature, 
and which ones are stable, we evaluate their total gravitational 
potential and thermal kinetic energies and compare their central 
densities with the theoretically predicted maximum stable cen- 
tral density of self-gravitating starless cores. 

In a first step, we calculated the gravitational potential en- 
ergy of the cores as 



GM 2 

-^grav — 71 ^vir 
K 



(12) 



with G being the gravitational constant, and using R = r 2 (Eq.|6] 
and Table[5]), M = M 2 the total gas mass enclosed in the N 2 
contour (Eq.|7J mass not explicitely listed here), and the profile 
shape correction factor Qf v i r = 0.75(±0.15) ( [Bertoldi & McKee] 
1992 ; Sipila e t al.|20li) . The thermal kinetic energy is calculated 
as follows: 
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Fig. 6. Stability of starless cores. Top: Ratio of gravitational po- 
tential to thermal kinetic energy vs. total gas mass (1.36 x M H ), 
both integrated within the N 2 column density contours (Table[5]). 
The lower dashed horizontal line marks the bounding limit of 
_ £grav/£therm = 1> while the upper line marks the state of virial- 
ization at -iL rav /£therm = 2.' Bottom: Estimated central density 
(see Sect. 5^) vs. total gas mass for the same starless cores. The 
dashed line marks the maximum stable density of a pressure- 
supported, self-gravitating modified (nonisothermal) BES (with 
photelectric h eating at the core bo undary taken into account) as 
calculated by |Keto & Casellij ( [2008] their Fig. 14). Typical esti- 
mated uncertainties on the source parameters are indicated as Y 
error bars on CB 4. The uncertainty in mass corresponds to the 
size of the label. 



with k being the Boltzmann constant, ji = 2.33 the mean molec- 
ular weight, mn the proton mass, and Af H x T being integrated 
within the N 2 contour. Figure[6](top) shows for the starless cores 
in our sample the ratio of gravitational to thermal energy vs. total 
gas mass M 2 . The turbulent energy content of the starless cores 
(derived from C 18 0(2-l) line widths; Lippok et al., in prep.) is 
by a factor 20 - 50 smaller than the thermal energy for all cores 
and is thus negligible in this balance. Given the uncertainties 
involved in this estimate, this energy ratio indicates that CB17 - 
SMM, CB26-SMM2, CB27, and B 68 are all gravitationally 
bound and approximately stable. CB244 - SMM2 must be su- 
percritical, which is consistent with the conclusion derived by 
|Stutz et aL] ( |2010 ), and CB 4 does not seem to be gravitationally 
bound. 

We also estimate the central volume number density of the 
starless cores by assuming spherical symmetry and dividing the 
peak column density No by the flat core radius r\ (Eqs.|3j [4] 
and Table[5]) and plot it in Fig.[6](bottom) against the total gas 
mass M 2 . Also indicated in this plot is the maximum stable cen- 
tral density of a pressure- supported, self-gravitating modified 



(nonisothermal) Bonnor-Ebert sphere ( BES;|Bonno r 1956, Ebert 
1955]), calculated by |Keto & CaselE] ( [200^ their Fig. 14). For 



comparison with the isolated globule cores, we adopt their re- 
sults for the model with photoelectric heating at the core bound- 
ary. This evaluation supports the previous conclusions even more 
clearly, i.e., CB17 - SMM, CB 26 - SMM2, CB 27, and B 68 are 
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all gravitationally bound but seem to be thermally subcritical and 
just at the limit of stability. 

This complies well with the notion that the mean radii of the 
globules are simi lar to or only slightly larger than their Jeans 
lengths (Sect 5.3). Interestingly, the projected separations be- 
tween the two sub-cores in CB 26 (0.15 pc) and CB 244 (0.09 pc) 
agrees quite well w ith th e typical Jeans lengths in these clouds 
(~ 0.13 pc; see Sect. 5.3), indicating that large globules tend to 
fragment gravitationally on the Jeans scale. Also interesting in 
this respect is that in both globules, one of the two cores has 
already formed a protostar, while the other one is still star- 
less/prestellar. 

CB4 is clearly a subcritical globule and must be purely 
pressure-confined. CB 244 - SMM2 must be a supercritical 
prestellar core at the verge of collapse, unless it has signifi- 
cant other nonthermal support like, e.g., strong magnetic fields. 
The CB 244 region has been studied by sub-mm polarization ob- 
servations reported in |Wolf et al.| ( [2003] ), and the SMM1 core 
has indeed been found to exhibit a relatively strong magnetic 
field (257 yuG), but the signal-to-noise for the SMM2 core turned 
out to be too low to derive meaningful field strengths. Hence, 
CB 244 - SMM2 is a good candidate to look for kinematic infall 
signatures or signs of a possibly already existing FHSC. 

Finding the majority of starless cores are bound and approx- 
imately stable is an expected result of a natural selection ef- 
fect. Due to the short collapse time scales (at most 10 4 yrs, e.g., 
Schonke & Tscharnuter 2011), unstable cores should be much 
rarer than stable cores which can have lifetimes of up to 10 6 yrs 
(e.g.,|Enoch et al.|2008|). 
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Fig. 8. Radially averaged column density profiles (bottom), 
LoS-averaged dust temperature (center), and luminosity-to-mass 
ratio profiles (top), for individual sources (thin black solid lines) 
as well as mean profiles (thick solid lines) of starless and proto- 
stellar cores. Note that the latter ones are fits to the mean profiles, 
while we list the means of the indivudual profile parameters in 
Table[5] which is not exactly the same. Dashed lines in the right- 
most panels indicate the 1 cr uncertainty ranges of the mean pro- 
files. 



5.6. Thermal structure of the cores 
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Fig. 7. Bolometric luminosity vs. total hydrogen mass of the en- 
tire globules (left), the envelopes only (without the dense c ores ; 
center) and the dense cores only (right; see Tableland Sect. |5.6| ). 
Filled quares in the right-most panel mark protostellar cores, 
empty circels mark starless cores, and the open square marks 
the core of CB 130 with the embedded VeLLO. 

Figure[7] shows the bolometric luminosity vs. total hydro- 
gen mass of the entire globules, the dense cores only, and 
the envelopes only (globule minus core). All globules turn 
out to have about the same luminosity-to-mass ratio of 1.1 ± 
0.6L o /M o , irrespective of whether they harbor starless cores or 
protostars. The mean luminosity-to-mass ratio of the envelopes 
(cores subracted) is very similar, with somewhat smaller scat- 
ter (1 ± O.5L /M ). A significant difference in the luminosity- 
to-mass ratios between star-less and protostar-harboring sources 
exists only for the central dense cores, with starless cores having 
0.5 + 0.25 L /M and protostellar cores having 2.5 + 1 L /M . 
This alone already points toward the dominating effect of the ex- 
ternally heated envelopes for the globules and the very local ef- 
fect of internal heating by embedded (low-mass) protostars. The 



embedded protostars have very little effect on the overall energy 
balance of their host globule. 

This dominating effect of external heating and shielding be- 
comes even more evident in the radial profiles of the dust tem- 
perature and the local luminosity-to-mass ratios. Figure[5] shows 
these radial profiles along with the column density profiles for 
all resolved starless and protostellar cores as well as the mean 
profiles for these two groups. Both starless globules and glob- 
ules harboring protostars approach at their outer boundaries, at 
radii of ~ 5 x 10 4 AU and column densities of ~ 10 21 cm -2 , dust 
temperatures of 15.5 + 1.5 K and local luminosity-to-mass ra- 
tios of 2 + 1 L /M . From the outer boundaries inward, both the 
dust temperatures and local luminosity-to-mass ratios decrease 
to reach mean values of 12- 13 K and ~ 0.4L o /M o at radii of 
~ 5000 AU. Hence, all of these cores are far from being isother- 
mal. They are heated externally by the IRSF and their inner- 
most dense centers are well-shielded from the energetic short- 
wavelength radiation of this heating source. Local heating by 
the embedded protostars raises the dust temperature and local 
luminosity-to-mass ratio only inside radii of ~ 5000 AU, where 
the derived local temperature is only a lower limit due to the 
effects of beam-convolution and LoS- averaging. However, there 
are also some systematic differences between certain globules, 
which we attribute to differences in the shielding by extended 
halos and difference in the local ISRF. This is being discussed 
in the next section. Table[5] also lists the column density and 
temperature profile fit parameters of the mean profiles of pro- 
tostellar and starless cores in our sample. Note that these mean 
profile parameters are averages of the individual source profile 
parameters, while the mean profiles shown in Fig. [8] are derived 
by fitting the averaged profiles. These observationally derived 
mean profile parameters can be used to construct representative 
cores for radiative transfer models. Such models must be able to 
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explain not only the outer and inner dust temperatures, but also 
the temperature contrast between core and envelope. 

5. 7. Constraints on the interstellar radiation field 
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Fig. 9. Central minimum (LoS-averaged) dust temperature, T- m , 
of starless cores vs. peak column density No. Representative er- 
ror bars are indicated only for CB 4. 

Figure[9] shows the central minimum (LoS-averaged) dust 
temperature, T- m , of all starless cores in our sample vs. their peak 
column density No (see Table[5]). The very tight anti-correlation 
between total column density and central dust temperature indi- 
cates very clearly that the inner temperatures of passive molec- 
ular cloud cores (without internal heating sources) are controled 
by the shielding against the ISRF. It also indicates that these 
nearby and isolated globules are all exposed to about the same 
ISRF. The correlated errors between T- m and No are much smaller 
than the range of values for these sources and therefore do not 
affect this conclusion. 

We also use the observed outer dust temperatures and halo 
column densities (N out ) of our sources to derive constraints on 
the mean ISRF. For this purpose, we compare the observed 
temperature profiles of the cores in this paper to a grid of ID 
radiative transfer models (Shirley et al., in prep.). Dust tem- 
peratu res are self-consi stently calculated using the CSDUST3 
code ( [Egan et aT]|1988| ) for an input density profile, input dust 
opacities, and input heating conditions. The density profiles are 
isothermal BES with the BES temperature derived iteratively 
at the half-central density point of the profile from the radia- 
tive transfer calculation. The central densities of the BESs vary 
from 9 x 10 3 cm 3 to 6 x 10 6 cm 3 with an outer radius set to 
30,000 AU. 

The dust opacities are taken from OH5 and are augmented 
by scattering coefficien ts according to the prescription given in 
Young & Evans | p005|). The heati ng is derived from the ISRF 
plotted in Figure 3 of |Shirley et al.| ([2005 ) and is modified in two 
ways: the overall strength of the UV to FIR ISRF is multiplied 
by a factor (SisrfX an d the model core is assumed to be embed- 
ded within a medium with extinction given by A v (R 0Ut ) affecting 
only the shortest wavelengths (A < few //m) of the ISRF. The 
final grid consists of 3600 models with self-consistently calcu- 
lated dust temperature profiles. 

Figure 10 compares the modeled relation between the outer 
temperature (at 30,000 AU), the relative strengths of the ISRF, 
and the additional shielding of the cloud against the ISRF by a 
halo for a model core that matches the mean column density pro- 
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Fig. 10. Observed mean outer dust temperatures of the glob- 
ules vs. column density of the outer halo (see Table[5]). The 
conversion to Ay was done with the mean A v -to-Nu ratio of 
(1.5 ± 0.5) x 10" 21 derived from the comparison of NIR extinc- 
tion maps (not used further in this paper) of 9 globules and the 
Nu maps derived in this paper and thus remains somewhat un- 
certain. Errorbars mark uncertainties of ±0.5 K for and 20% 
for Af H (out), repectively. Plotted as grayscale and contours is 
the dust temperature at the outer radius (30,000 AU) of a model 
cloud, derived with the radiative transfer modeling descibed in 
Sect. |5.7| versus shielding of the cloud against the ISRF by a 
halo and as function of the relative strengths of the ISRF. 



file of the observed starless globules (Fig. [8]) with the observed 
outer dust temperature and column densities of the outer ha- 
los (see Table[5]). The observed outer temperatures range from 
13. IK to 19. 5 K and occupy the upper half of the model outer 
temperatures which range from 8.1 K to 17.8 K. The general 
shape of the observed dependency of the outer temperature on 
the shielding against the ISRF agrees well with the model calcu- 
lations, suggesting that the outer temperature is indeed controled 
to first order by the shielding against an approximately uniform 
(not necessarily isotropic) ISRF. The only source that deviates 
significantly from the observed relation toward higher outer tem- 
perature or stronger ISRF is BHR 12 (CG30). This cometary 
globule in the Gum nebula region is exposed to the Vela pulsar 
and supernova remnant (Larg e et al.|19 68 ) and to the 04 super- 
giant star f Puppis, and might thus inde ed ex perience a stronger 
local ISRF than the other globules (Fig. |A.6| ). 

However, the comparison between observations and model 
calculations also shows that the isolated globules hav e on aver- 
age highe r outer dust temperatures than predicted by Shirl ey et] 
[aT] ( |2005| ). This is probably due to a combination of a stronger 
ISRF than adopted and lower submm opacities for grains with a 
smaller size distribution than OH5. For this calculation, we used 
the same dust opacity model for both the dense cold cores and 
for the thin extended envelopes. OH5 opacities, which are calcu- 
lated for coagulated dust grains with ice mantles, are certainly 
not appropriate for the thin extended halos around the glob- 
ules. Though opacity models for smaller grains will not signif- 
icantly affect the observational derivation of dust temperatures 
from the SEDs, they would lead to smaller values of the derived 
outer column densities (by a factor of a few). At the same time, 
smaller grains would lead to higher outer dust temperatures in 



17 



R. Launhardt et al.: The structure of low-mass cloud cores 



the self-consistently calculated models. These issues require fur- 
ther modeling before robust conclusions on the strength of the 
ISRF can be drawn and will be addressed in a forthcoming pa- 
per that will also include information from NIR extinction maps. 

Nevertheless, the observed relation between outer dust tem- 
peratures and halo column densities complies very well with the 
general picture of external heating by the ISRF and, at the same 
time, suggests that our backgro und fl ux level treatment and ex- 
tended emission recovery (Sect. |4.1| ) worked well and that the 
detected "halos" indeed belong to the globules and do not repre- 
sent unrelated background material. 

5.8. Environmental effects 

Although all globules selected for this study appear at first 
glance compact and not very filamentary, eight out of the 12 
globules are cometary-shaped with one sharp rim and either a 
long cometary tail seen both in the optical cloudshine as well as 
in the FIR dust emission (e.g. CB 6, CB 17, CB 68) or at least a 
very diffuse side opp osite the shar p rim (e.g., CB 26, CB 230; 
see Tableland Figs. A.l through [A. 12 ). Most of these glob- 
ules where not classified originally as cometary globules like, 
e.g., BHR12 (CG30; |Zealey et al][T983] ), but appear at first 
glance rather round and isolated when seen on less deep opti- 
cal images or older dust emission maps. The diffuse tails have 
projected physica l leng ths between 0.3 pc (CB 27; Fig. |A.5| ) and 
2.6 pc (CB 6; Fig. A.2 ). In some cases, like CB 68, the long dif- 
fuse tail connects the globule to larger (diffuse) clouds (Figs. [2] 
and |A.7| ). In other cases, like BHR 12, the cometary tail is point- 
ing away from an energetic irradiating source (Fig. A.6 ). In some 
sources, like, e.g., CB 17, neither a connection of the long tail 
to another structure nor an illuminating source could be clearly 
identified. This diversity suggests that cometary shape and tails 
may not have the same physical origin in all sources. 

We have also evaluated the location and orientation of warm 
rims on the outer edges of the globules (like, e.g., CB 26: Fig. [27]) 
and their relation to (potential) external heating sources, like 
nearby stars or the galactic plane. Even though the one-sided 
warm rim appears at first glance questionable in some cases 
because we have recorded no information for the temperature 
of the opposite rim (e.g., CB4 and CB 17: Figs.[24l and |26|), 
the existence and location of the one-sided warm rims is in- 
deed significant. The warm rims, were present, are on average 
2 ± 0.4 K warmer than the opposite cold rims. Their presence can 
be understood by inspecting the PACS and SPIRE maps of, e.g., 
CB 4 (Fig. 12). While the long-wavelength SPIRE maps, which 
are to first order sensitive to column density, show an approxi- 
mate symmetry between northern and southern side of the glob- 
ule (apart from the cometary tail), the shorter- wavelength PACS 
maps, which are more sensitive to temperature, show a signif- 
icant asymmetry with much stronger emission on the southern 
side. Indeed, after fitting the data, the column density contours 
indicate an approximate symmetry between northern and south- 
ern side down to a few times 10 20 cm -2 (Fig. 24 ). At lower col- 
umn densities, we can still derive the temperature at the warm 
southern rim where the short- wavelength PACS emission is suf- 
ficiently strong. But on the northern side, the short-wavelength 
PACS flux densities drop below the detection threshold. Hence, 
we can no longer reliably derive a temperature for the cold ma- 
terial which is still traced by the SPIRE maps and may have the 
same column density as the warmer material on the southern 
side. When allowing the fitting routine to derive temperatures 
from only three (SPIRE) wavelengths (cf. Sect. 4.2 ), we obtain 
a temperature map that fully confirms the above result and rules 



out a misinterpretation due to a detection bias driven the data re- 
duction. However, these SPIRE-only derived temperature maps 
are too unreliable and noisy to be used in the analysis of the 
source structure presented in this paper. Furthermore, the com- 
parison of the results for CB 244 obtained with different ver- 
sions of the data reduction and calibration pipeline (Sect. |6.1| ) 
also show that the extent and location of the warm south-eastern 



rim in CB 244 (Fig. 35 ) remains unaffected by the specific data 
reduction version or background subtraction level. 



Figures A.l through A.12 show the large-scale morphology 
and surroundings of all regions studied in this paper, including 
the schematic indication of the warm rims and the location of 
potential external illumination sources. We find no clear indica- 
tion for a general link between the loaction of the warm rim side 
in a globule and the location of identified potential external heat- 
ing sources. In some globules, the warm rim points toward the 



galactic plane ( CB 26, B 6 8, CB 68, CB 230, CB 244: Figs.|A4 
A.8||A.7||A.11| and |A.12| ), in others it points in the opposite di- 



rection (CB 4, CB 27, CB 130: Figs.[AT][A3) and|A^9]). In CB 6 
(Fig. |A.2| ) and CB 17 (Fig. |A.3| ), there are bright nearby stars 
identified that could be located at the same distance and thus 
illuminate the globules directly. However, the warm rims do not 
point toward the potential illuminators. Except for BHR 12, no 
other nearby heating or illuminating sources could be uniquely 
identified. However, even in this globule, the tail pointing away 
from the illuminating source is significantly warmer than the 
side facing it (Figs.[29| and |A.6| ). Hence, we have to conclude 
that, as with the cometary tails, the link between warm rims and 
illuminating sources or the physical reason for the asymmetric 
warm rims remains inconclusive and needs to be addressed by 
future studies. 



6. Discussion 

6. 1. Calibration and data reduction uncertainties 

The total flux sc ale u ncertainty of the Herschel data is assumed 
to be 15% (Sect. |4.1| ), while that for the ground-based submm 
bolometer data is about 20% (Sect. |3.2| ). These flux calibra- 
tion uncertainties propagate directly into the column density and 
mass estimates. However, the effect on the dust temperature es- 
timates should be much smaller for the following reason. Since 
all bands within one of the Herschel instruments were photo- 
metrically calibrated in a consistent way, any systematic uncer- 
tainties in the effective temperature of the calibrators should be 
correlated between the bands and thus to first order not affect 
the flux ratios that determine the temperature. There could still 
be a systematic offset between PACS and SPIRE due to differ- 
ent calibration procedures. However, if such an offset exists, it 
would have become immediately evident in the combined pixel- 
by-pixels SEDs since we have more than one wavelength band 
from each instr ument. Since, after the background offset subtrac- 
tion (Sect. |4.1| ), we have not noticed any such offset, it must be 
very small (see also the integrated SEDs in Fig. [4]). Similarly, 
any systematic errors that we could have possibly i ntrod uced 
in the various post-processing steps described in Sect. |4.1| (e.g., 
spatial filtering) should also be correlated to first order between 
the Herschel bands. Here we assess these effects quantitatively 
and estimate the uncertainties introduced by using different and 
older versions of HIPE and Scanamorphos for different sources 
(Sect. |3.1| ). I particular, we have reprocessed all Herschel data 
from one source (CB 244, originally reduced with HIPE v. 5/6 
and Scanamorphos v. 9) with the newest versions of HIPE (v. 9) 
and Scanamorphos (v. 18), with all additional post-processing 
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steps (e.g., background subtraction) unchanged. The new reduc- 
tion resulted in up to 20% lower PACS fluxes (as compared to 
the old reduction), but left the SPIRE fluxes unchanged within 
« 1%. As a result, the LoS -averaged dust temperatures derived 
with the newly processed Herschel maps are on average 2.2% 
lower than the old values, and the resulting column densities are 
on average 7.5% higher. We also interpret these numbers here as 
a general assessment of the possible changes that new data re- 
duction versions might bring, since HIPE v. 9.1.0 may not be the 
"final" version. 

The complementary ground-based submm data were ob- 
tained independently with different instruments and processed 
with the respective instrument pipelines. Not only are the cali- 
bration uncertainties between the instruments and with respect 
to Herschel no longer correlated, these data are also much more 
and in a different way affected by spatial filtering. The sum of 
these effects would make any results based on simple submm 
flux ratio maps very susceptible to a number of calibration un- 
certainties. However, we use all data in a combined analysis in 
which the dust temperature estimate is to first order determined 
by the flux ratios of the Herschel bands, while the column den- 
sity is to first order determined by the long-wavelength SPIRE 
bands and the submm data. In order to quantify the effect of the 
calibration uncertainties on our temperature estimates, we have 
carried out a Monte Carlo test on one pixel in the center of B 68, 
artificially increasing all Herschel fluxes by 15% and randomly 
modifying the three individual submm fluxes with a rms of 20% 
and fitting the SEDs in the same way as we did for the normal 
analysis. As a result, the mean temperature of that pixel was ele- 
vated by 0.2 K from 1 1 .8 K to 12.0 K with a rms scatter of 0.2 K. 
The mean column density was elevated by 6% with a rms scatter 
of 10%. 



The background subtraction described in Sect. |4.2| could 
have potentially introduced additional uncertainties. However, 
since 4g ^ 2? V (10K) at all wavelengths for all our sources, 
thanks to the initial source selection in regions of low back- 
ground, the difference in the final temperature estimates between 
considering 4g or setting 4g = is AT" < 0.15 K. The uncer- 
tainty introduced by not knowing the exact level of the back- 
ground, but still subtracting a mean guess of 4g is correspond- 
ingly smaller, i.e., negligible. We therefore estimate the over- 
all relative 1 <x uncertainty introduced by flux calibration and 
data processing uncertainties to be cr T « 3% and cr N « 15% 
(quadratic sums of individual contributions described above). 

6.2. Uncertainties and limitations of the modeling 
approach 



Througho ut this paper, we used OH5 opacities (Qssenko pf & 
Henning 1994j ) for all sources and regions. Although this opac- 
ity model seems to be a reasonable assumption for this kind of 
sources, it is not well-constrained observationally and may not 
accurately reflect the actual dust properties of the individual ob- 
jects. Furthermore, a single opacity model may not apply to both 
the dust in the cold dense cores, where coagulation and ice man- 
tle growth may have affected the effective optical properties of 
the grains, and that in the thin and warm outer envelopes. Since 
the data do not provide a direct and unambiguous constraint on 
the dust opacities (this will be subject of a subsequent paper), we 
tested the dependency of the temperature and column density es- 
timation on different opacity models and report the results here 
in the form of an uncertainty discussion. In particular, we test 
for one source (B 68) how the temperature and column density 
estimates derived with the following five dust opacity models 




R [arcsec] 

Fig. 11. Comparison of radial column density profiles of B 68, 
derived with the LoS-averaged temperature recovery method de- 
scribed in this paper (black solid line, uncertainty range marked 
by dotted lines), with the assumption of isothermal dust emission 
(dash-dotted line), and with a ray-tracing approach that takes 
LoS temperature gradients into account (gray dashed line, uncer- 
tainty range marked by gray- shaded area, Niel bock et al.|2012 
corrected for OH5 opacities). 



differ with respect to the use of OH5 opacities: (1) OH1 (MRN 
size distribution, Mathis et al.|1977] ), (2) OH8 (coagulated grains 
like OH5, but with thick ice mantles), (3) Weingartner & Draine 
opac ities for R v = 5.5 (Weingartner & Draine 2001 Draine 
|2003p ] and (4) two pure beta models with £ = 1.5 and£ = 2.5, 
normalized to OH5 (which has ft « 1.9) at A = 500 fim. The 
first three of these models lead to deviations in the derived tem- 
perature w.r.t. OH5 of about 4% (0.4 K at 10 K), with the maxi- 
mum deviation occuring for WD01 in the warm envelope (9% or 
1.6 K lower than the 18 K derived with OH5). The two extreme 
beta models result in « 10% higher temperatures (w.r.t. OH5) for 
J3 = 1.5 or 14-18% lower temperatures for/3 = 2.5, respectively. 
Hence, we can safely assume that the systematic uncertainty of 
the dust temperatures listed in Table[5] introduced by the uncer- 
tainty in the dust opacity model is AT « ± 1 K. Since the as- 
sumption of ISM-like dust (e.g., OH1), which is probably more 
appropriate for the envelopes than OH5, leads to a reduction of 
the derived temperature by only a few per cent (w.r.t. OH5), 
the LOS-averaged temperature contrast between cores and en- 
velopes could be systematically overestimated by « 0.5 K, i.e., 
less than the absolute temperature uncertainty. 

The LoS-averaged dust temperature maps presented in this 
paper provide a robust estimate of the actual dust temperature 
only for the envelopes in the projected outer regions, where the 
emission is optically thin at all wavelengths and LoS temper- 
ature gradients are negligible. Toward the core centers, where 
cooling and shielding or embedded heating sources can produce 
significant LoS temperature gradients and the observed SEDs 
are therefore broader than single-temperature SEDs, the local 
central dust temperature is overestimated (in the case of a pos- 
itive gradient in cold sources) or underestimated (negative gra- 
dient in internally heated sources). E.g., for B 68, where we al- 
ready derived a model that accounts for LoS gra dients, we de- 
rive a local central dust temperature of 8.2^ K (Nielbock et al. 
2012), which is 3.7 ± 1 K lower than the LoS-averaged temper- 
ature of 1 1.9 + 1 K we derive in this survey overview paper. For 
CB 17, where the outer dust temperature is lower than in B 68 
(Table|5]), we derive a local central dust temperature that is only 



See http://www.astro.princeton.edu/~draine/dust/dustmix.html 
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« 1 K lower than the LoS -averaged dust temperature maps de- 
rived in this paper (Schm alzl et al.||20T3 ). Dust opacity model 
uncertainties cancel out to first order if we compare only results 
obtained with the same opacity curves. I.e., in this paper we sys- 
tematically overestimate the minimal inner dust temperatures by 
1 - 4 K, underestimate the outer envelope dust temperature by 
< 0.5 K, and underestimate the temperature contrast between 
cold core and envelope by 1 - 4 K. 

The effect of systematically overestimating the actual min- 
imum inner dust temperatures on the column density and mass 
estimates is also not negligible and can be quite significant in 
individual cases. Figure TT| compares the column density pro- 
files of B 68 derived with an isothermal assumption, the LoS- 
averaged SED fitting used in this paper, and the full RT approach 
used in Nielbo cF et al. | pOT2| ) . In this rather extreme case (one 
of the highest core - envelope temperature contrasts in our sam- 
ple), the LoS-averaged SED fitting underestimates the peak col- 
umn density by a factor of « 1.7 and the total mass by « 15%. 
In sources with less-pronounc ed temperature contra st between 
core and envelope, like CB 17 ( Schm alzl et al.|2013] ), these sys- 



tematic effects are smaller. Still7 Fig.| 1 1 1 also illustrates that the 
column density profiles recovered by our SED fitting approach 
are a much better approximation to the actual underlying source 
structure than those derived with an isothermal assumption. 



The local temperature minimum and column density peak 
values are also affected by the angular resolution of our anal- 
ysis. The 36'.' 4 beam size we use for all data may smooth out 
the central parts of the more centrally peaked profiles in some 
sources. Although we do not consider this effect being very sig- 
nificant for the starless cores, since the derived flat c ore radii are 
all significantly larger than the beam (see Sect. |5.3| ), it can only 
be verified quantitatively by using all data in their original res- 
olution and employing full forward modeling. Beam- smoothing 
does affect the temperature and column density estimates much 
more at the position of embedded protostars where temperature 
and density gradients are far more pronounced and definitely not 
resolved. Here, the local temperatures at the position of embed- 
ded protostars are usually strongly underestimated and the tem- 
perature and column density maps are often corrupted by even 
very small residual pointing mismatches. 

Hence we can summarize that the uncertainties introduced 
by calibration and data reduction issues amount to « 3% for the 
temperature and « 15% for the column density. The uncertain- 
ties related to the dust opacity models are significantly larger and 
amount to ±1 K (« 10%) for the temperature and a factor of up 
to 3 for the column density and mass. Ignoring LoS temperature 
gradients leads to an overestimation of the central temperature 
minima in starless cores of « 1 - 3 K and an underestimation of 
the central column density by up to 50%. The temperatures and 
column densities of the outer envelopes are much less affected 
by this LoS averaging. 

The largest (systematic) uncertainties of the derived temper- 
atures and (column) densities are clearly related to our uncertain 
knowledge of the underlying dust mass absorption coefficient 
and its spectral dependence. Therefore, our goal in future work 
will be to exploit the Herschel data presented here along with 
additional complementary data to derive specific observational 
constraints on the properties of the dust in these globules (see 
also end of Sect. [7]). 



6.3. Comparison with earlier observations and models of 
dense cores 

In this section we discuss some examples of previous work based 
on resolved observations of starless low-mass cores. For com- 
parison with our new Herschel results, we specifically discuss 
published measurements of density and temperature profiles, as 
well as different attempts to model said temperature and density 
profiles. 

The observed nonisothermality of dense clouds, derived 
from Herschel data, was already reported by other authors (e.g., 
Stutz et al.||2010| |Peretto et al.|| 2010| JBattersby et al.|[20TT| 

|20l3) . While the |Stutz et al.| 



Wilcock et al.||2012l 

JL ' 



IPitann et al 



( 2010 ) p aper was the pilo t study for the globule survey presented 
here,|Peretto et al.|(|2010j),|Battersby et al.| ( |201l] ), | Wilcock et al.| 
( |2012| ), and |Pitann et al.| ( |2013| ) found that the LoS-averaged dust 
temperatures in IRD Cs also systematically drop toward the high- 
est column densities. Wilcock ^et al.| ( |20l21 Fig. 5) found temper- 
ature profiles of IRDCs that are qualitatively very similar to the 
T profiles of the starless globules we report in this paper (Fig. [3]), 
although their outer temperatures were on average ~ 5 K higher 
than those found here (20 K vs. 15 K). This is not surprising be- 
cause IRDCs are generally located in distant regions of high- 
mass star formation with strong MIR background and may on 
average be exposed to a stronger ISRF than the nearby isolated 
globules in the local spiral arm. In more nearby sources, low cen- 
tral temperatures and positive temperature gradients were also 
inferred from observations and m odelin g of molecular lines in 
prestellar cores (e.g., Paga ni et al.|2007 Crapsi et al.|2007 ). 



Ward-Thompson et al.| (|2002|) already found observational 



evidence for temperature gradients and colder interriors of star- 
less cores using flux ratio maps of 17 and 200 yum observed 
with ISOPHOT. |Bacmann efaTl (|2000) observed 24 dense star- 
less cores in absorption with ISOCAM at 6" angular resolution. 
From modeling of the 7 fim absorption against the MIR back- 
ground, they derived core structures characterized by a central 
flat region of 7?fl at « 4000 - 8000 AU, central column densi- 
ties of a few 10 22 cm -2 , and o uter radial column den sity profiles 
close to Nu oc r" 1 . While the Bac mann et al.| (|2b00) mean peak 
column densities agree well with the values we derive, our anal- 
ysis results in somewhat larger flat core radii and significantly 
steeper column density profiles (Table|5]), which might at least 
partially be attributable to the larger beam sizes of our maps as 
compared to the Bacmann work. 



Schnee & Goodman (2005) used 450 and 850 ^m SCUBA 



observations to derive the dust temperature profile of the nearby 
starless core TMC-1C from the 450/850 /mi flux ratio (assum- 
ing optically thin emission and a fixed ft = 1.5). They found a 
positive temperature gradient with a central dust temperature of 
6 K and an outer temperature of 12 K at the edge of the core (at 
0.08 pc). They derived similar column densities as we find for 
the globule sources and modeled the density profile with a bro- 
ken power-law with a shallow density profile inside ~ 4000 AU 
and oc r -1 - 8 outside. Since the break radius of a broken power- 
law model of the density is not the same as the parameter r\ of 
our column density profile fits (Eq.|3]) and should be somewhat 
smaller, their flat core radius is comparable or slightly smaller 
(up to a factor of two) than what we derive for the globule 
cores. Their outer density profile exponent of 1.8 is significantly 
smaller than what would be implied by the column density pro- 
files we derive here (Table|5]) and the density profile of B 68 we 
derive in |Nielbock et al.|(2012| ). However, both the dependence 
of the resulting temperature on the above-mentioned assump- 
tions as well as the non-negligible calibration uncertainties of 
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the SCUBA data, in particular at 450 fim, make the temperature 
scale and the (column) density profile derived from such flux 
ratios rather uncertain. 



|Shirley et al] ( [2QQQ| ) also observed with SCUBA the 450 and 
850 //m dust emission of 21 low-mass cores in different evolu- 
tionary stages. Based on these observations, |Evans et al. ( 2001| ) 
modeled the emission of three starless cores within the context 
of a self-consistent radiative transfer analysis assuming a non- 
isothermal BES, modified OH5 dust opacities, and a modified 



Black| ( ^994| ) and |praine| ( [T978| ISRE They found central den- 
sities between 10 D and 10° cm" 3 and required only fractions 
(0.3-0.6) of the adopted ISRF to reproduce the 450/850 /im 
flux ratios. They also derived positive temperature gradients with 
central dust temperatures of 7 - 8 K and outer temperatures of 
12 - 13 K and found that the ISRF was the dominan t heating 
source in the absence of a nearby star (see also Zucconi et al. 
2001). Given the fact that our LoS -averaged dust temperatures 



overestimate the actual central temperatures by 1 - 4 K, their de- 
rived central temperatures are still 2 - 3 K lower than what we 
find here. However, one has to keep in mind that this is based on 
very low number statistics and that their three cores are on aver- 
age denser and more massive than ours. Therefore, this small dif- 
ference in central temperatures might be real and reflect different 
source properties. Their outer temperatures, which are not well- 
constrained by the observations, are also 2-3 K below the obser- 
vationally well-constrained values we derive. However, without 
comparing the same sources, it is not possible to verify whether 
this reflects real differences in the shielding or is a result of the 
specific assumptions for the dust properties and the strength of 
the ISRF they made (see also discussion in Sect. 5/7). However, 
the same remarks apply as made above about the uncertainty of 
such quantities that are based only on SCUBA flux ratios. 

|Shirley et al.| p002| ) modeled the Class cores from the 
|Shirley et al.| ( |2000| ) sample. In order to reproduce the 450 and 
850/mi fluxes, they included the internal radiation field from 
the protostar in addition to the heating from the ISRF. They 
found that the dust temperature exceeded 20 K only in the in- 
ner * 10 3 AU, reached * 10 K at r * 5000 - 8000 AU, and 
increased again slightly toward the outer boundary of the core at 
a few 10 4 AU. This is consistent with what we find, including the 
temperature rise at r > 10 4 AU, taking again into account that we 
overestimate the inner minimum temperature by up to 3 K, and 
that their models appear to underestimate the outer temperatures 
by 2 - 3 K. 



Stamatellos et al. (2005 ) calculated temperature profiles in 



the context of an evolutionary model, starting with a 5.4 M© BES 
and assuming OH5 dust opacities and a modified Black (1994) 
ISRF. In the initial stage, when only the ISRF heats the core, they 
found a ce ntral temperature of ~7 K and outer temperatures of up 
to -20 K. |Stamatellos et aT]p007| ) further studied the influence 
of external heating and shielding on the temperature profile of 
cores in the p Oph cloud, taking into account also the surround- 
ing molecular cloud extinction. They found again that the inner 
core temperatures (at r < 10 4 AU) were always below 10 K, but 
that they reached rim temperatures (at r « 2 x 10 4 AU) of only 
-13-15 K. They noted that in Taurus, where the ambient radia- 
tion field is weaker, the core temperatures are found to be some- 
what lower with a central value of « 7 K and rim temperature of 
« 12 K for TMC-1C. These temperatures and temperature con- 
trasts are again consistent with what we now measure from the 
Herschel data in similar cores. In particular the outer tempera- 
tures we measure for the isolated globules with their thin halos 
are all inbetween the values they derived for the unshielded and 
the well-shielded cores. 



7. Summary and conclusions 

We have studied the thermal dust emission from 12 nearby iso- 
lated globules using FIR dust continuum maps at five wavelength 
bands between 100 and 500 jim obtained with the Herschel satel- 
lite and submm continuum maps at up to three wavelengths 
beween 450 /mi and 1.2 mm obtained with ground-based tele- 
scopes. Great care was taken in both the selection of truly iso- 
lated and previously well-characterized sources in regions of ex- 
ceptionally low background emission and in restoring and cali- 
brating all maps in a consistent manner to enable a combined in- 
depth analysis of all data. This careful preselection, data treat- 
ment, and the very deep 100 //m maps (~ 50/Jy/arcsec 2 rms) 
turned out to be essential to enable our analysis method and lead 
to the results presented here. 

We calibrate all data to the same flux scale, subtract all gen- 
eral background levels in a constistent way, convolve all maps 
to the same angular resolution, assign weights to each pixel ac- 
cording to calibration uncertainties and rms noise, and extract 
full SEDs from 100 yum to 1.2 mm for each image pixel. These 
individual SEDs are independently fit (least squares) with single- 
temperature modified blackbody curves, using OH5 opacities for 
mildy coagulated dust grains with thin ice mantles. Thus we ob- 
tain for each image pixel a column density and an LoS -averaged 
dust temperature. The resulting dust temperature maps have rms 
noise levels of only a few tenth of a K and systematic uncer- 
tainties of less than + 1 K. We then analyze the dust tempera- 
ture and column density maps along with the integrated SEDs 
of all sources, including additional archival data at shorter wave- 
lengths, and obtain the following main results: 

1 . The globules in our sample have a mean distance from the 
Sun of 240 ± 100 pc, mean total gas mass of 7.3 ± 6.8 M , 
mean bolometric luminosity of « 6.7 L , and mean outer di- 
ameter of 0.45 ± 0.2 pc. Typical peak column densities (aver- 
aged over the 36'.'4 beam) are (Nu) ~ 3 x 10 22 cm" 2 , which 
corresponds to peak visual extinctions of about 45 mag. 

2. The dense cores of the globules have a mean total gas mass 
of 1 .6 ± 1 .2 M© and a mean diameter of 0. 12 ± 0.05 pc. They 
constitute on average 25 ± 6 % of the total globule mass. 

3. Out of the 12 globules in our sample, four contain one star- 
less core each, six contain a protostellar core (Class or I), 
and two globules contain each one starless and one protostel- 
lar core. 

4. We confirm the earlier classifications of protostars and em- 
bedded YSOs, but the better- sampled thermal SEDs do not 
allow for better evolutionary stage determination from the 
SEDs alone since projection effects dominate the spectral ap- 
pearance of even Class protostars. However, the Herschel 
data allow us now to also include starless cores in the 
L smm / Lboi vs. Tboi diagram (Fig. [5]). We find that all star- 
less cores have r bo i < 25 K and 10% < L smm / L bo i < 30%. 
In particular the criterion 7b i < 25 K can be used to identify 
isolated and externally heated starless cores, since even an 
embedded VeLLO will elevate r bo i above 30 K. 

5. Most globules are surrounded by a thin (A v ~ 1 - 3 mag) 
and warm (T d « 15-20 K) outer halo, the full extent of which 
we could not properly recover in most cases. These halos 
are evident in the dust emission-based column density maps, 
NIR extinction maps, as well as in many cases on optical 
images as cloudshine. 

6. The thermal structure of all globules is dominated by ex- 
ternal heating through the ISRF, with warm outer envelopes 
(15.5 + 1 .8 K) and colder interiors (< 13 + 1 K). We find a very 
tight anti-correlation between the minimal central tempera- 



21 



R. Launhardt et al.: The structure of low-mass cloud cores 



tures in starless cores and the total column density toward the 
outer edge that provides shielding against the ISRF (Fig. [9]). 

7. The outer dust temperature of the globules is to first order 
determined by external heating from an approximately uni- 
form ISRF and shielding by the thin extended halos. We find 
a clear anti-correlation between the outer dust temperatures 
and the mean column density / Ay of the halos (Fig. 10). 
Only one of the globules is exposed to a significantly higher 
than average local ISRF (BHR 12). 

8. The embedded (low-mass and low-luminosity) protostars 
raise the local temperature of the dense cores only within 
radii out to about 5000 AU, but do not significantly affect the 
overall thermal balance of the globules. 

9. Our LoS-averaged dust temperatures are accurate to about 
± 1 K at the outer boundaries of the globules (accounting 
for both statistical and systematic errors), but systematically 
overestimate the local central temperature minima by 1 - 4 K, 
depending on the temperature contrast to the outer rim. I.e., 
these globules have central temperatures in the range 8 - 12 K 
and temperature contrasts to the outer rim of 3 - 8 K. 

10. Out of the six starless cores in our sample, five cores are 
found to be thermally stabilized against gravitational col- 
lapse. Of these, four cores (CB 17-SMM, CB 26 - SMM2, 
CB 27, and B 68) have central densities that are very close 
to the maximum stable density of a pressure- supported, self- 
gravitating cloud, i.e., we cannot distinguish wheter they are 
starless stable or prestellar. CB 4 is thermally sub-critical and 
must be a purely pressure-confined starless core. The starless 
core in CB 244 - SMM2 is clearly super-critical and should 
very likely be a prestellar core at the verge of collapse. 

1 1 . The radius of the single-core globules is very similar to or 
only slightly larger than their Jeans length (~ 0.13 pc). The 
projected core separation in the two double-core globules 
agrees well with the Jeans lengths in these clouds, indicat- 
ing that large globules tend to fragment gravitationally on 
the Jeans scale. 

12. Many of the globules have one side that is signifcantly 
warmer at the outer rim compared to the other side. Although 
in some cases the warm side faces the direction toward the 
Galactic plane, we could not identify a consistent relation be- 
tween warm side and potential external illuminating sources. 



As a next step, we are currently analyzing individual starless 
cores by employing a ray-tracing algorithm that is able to restore 
LoS temperature gradients and derive the actual volume desnity 
structure of the cores by starting from a 1-D mode l and allow- 
ing for local deviations from spherical symmetry (Nielbo cTTet| 
al.|2012[|Schmalzl et al.|2013| ). This way we hope to better and 



most robustly constrain the dust temperature structure of these 
cores without introducing too many model-dependent assump- 
tions into the analysis. We are also working toward better con- 
straining the dust opacity models. For this purpose, we have al- 
ready obtained deep NIR extinction maps that recover the same 
angular resolution as the Herschel data and are currently ob- 
taining additional deep submm and mm continuum maps. We 
are planning to employ and compare different approaches, in- 
cluding SED fitting and ray-tracing like in this paper and in 
Nielbock et al.|(|2012|), self-co nsistent radiative transfer model- 



ing like in 



like in Kel 



Shirle y et al. 
y etal.|p012 



(2011 



. Furt] 



), as well as Bayesian approaches 
thermore, we are currently also us- 
ing complementary molecular line data to derive constraints on 
the distribution of turbulence and molecular gas-phase freeze- 
out onto dust grains in the sources studied in this paper (Lippok 
et al., in prep.). 
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Fig. 12. CB 4 (starless core): Visual (red) DSS2 image and Herschel FIR maps (color; log scale) at 100, 160, 250, 350, and 500 yum, 
with 1.2 mm continuum contours overlaid (20 and 40mJy/30 // beam). Herschel beam sizes are indicated as gray circles in the lower 
right corners. 
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Fig. 13. CB 6 (Class I YSO): Visual (red) DSS2 image and Herschel FIR maps (color; log scale) at 100, 160, 250, 350, and 500 /mi, 
with 850 yum continuum contours overlaid (180 and 300 mJy/20" beam). Herschel beam sizes are indicated as gray circles in the 
lower right corners. 
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Fig. 14. CB 17 (starless core and embedded Class I YSO): Visual (red) DSS2 image and Herschel FIR maps (color; log scale) at 
100, 160, 250, 350, and 500 yum, with 1.2 mm continuum contours overlaid (20, 40, and 65 mJy/20"beam). Herschel beam sizes are 
indicated as gray circles in the lower right corners. 
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Fig. 15. CB 26 (starless core and embedded Class I YSO): Visual (red) DSS2 image and Herschel FIR maps (color; log scale) at 
100, 160, 250, 350, and 500 yum, with 1.2 mm continuum contours overlaid (10, 26, and 150mJy/30"beam). Herschel beam sizes 
are indicated as gray circles in the lower right corners. 
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Fig. 16. CB 27 (starless core): Visual (red) DSS1 image and Herschel FIR maps (color; log scale) at 100, 160, 250, 350, and 500 yum, 
with 1.2 mm continuum contours overlaid (40 and 70mJy/30"beam). Herschel beam sizes are indicated as gray circles in the lower 
right corners. 
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Fig. 17. BHR 12 (Class 0/1 protostars): Visual (red) DSS2 image and Herschel FIR maps (color; log scale) at 100, 160, 250, 350, and 
500 yum, with 850/mi continuum contours overlaid (SCUBA; 120, 300, and 900mJy/15 // beam). Herschel beam sizes are indicated 
as gray circles in the lower right corners. 
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Fig. 18. CB 68 (Class protostar): Visual (red) DSS2 image and Herschel FIR maps (color; log scale) at 100, 160, 250, 350, and 
500 yum, with 850yum continuum contours overlaid (150, 300, and 600mJy/20 // beam). Herschel beam sizes are indicated as gray 
circles in the lower right corners. 
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Fig. 19. B 68 (starless core): Visual (red) DSS2 image and Herschel FIR maps (color; log scale) at 100, 160, 250, 350, and 500 yum, 
with 870 yum continuum contours overlaid (40, 120, and 220mJy/27 // beam). Herschel beam sizes are indicated as gray circles in the 
lower right corners. 
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Fig. 20. CB 130 (Class protostar): Visual (red) DSS2 image and Herschel FIR maps (color; log scale) at 100, 160, 250, 350, 
and 500 yum, with 1.2 mm continuum contours overlaid (20, 60, and 100 mJy/20" beam). Herschel beam sizes are indicated as gray 
circles in the lower right corners. 
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Fig. 21. B 335 (Class protostar): Visual (red) DSS2 image and Herschel FIR maps (color; log scale) at 100, 160, 250, 350, and 
500 yum, with 850 yum continuum contours overlaid (100, 300, and 900mJy/20 // beam). Herschel beam sizes are indicated as gray 
circles in the lower right corners. 
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Fig. 22. CB230 (Class I YSO): Visual (red) DSS2 image and Herschel FIR maps (color; log scale) at 100, 160, 250, 350, and 
500 yum, with 850 /mi continuum contours overlaid (70, 200, and 700mJy/20"beam). Herschel beam sizes are indicated as gray 
circles in the lower right corners. 
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Fig. 23. CB 244 (Class protostar and starless core): Visual (red) DSS2 image and Herschel FIR maps (color; log scale) at 100, 
160, 250, 350, and 500 yum, with 1.2 mm continuum contours overlaid (30, 80, and 180mJy/20"beam). Herschel beam sizes are 
indicated as gray circles in the lower right corners. 
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Fig. 24. Dust temperature map of CB 4 (color scale), overlaid 
with contours of the hydrogen column density (white: 10 21 cm" 2 , 
black: 0.2, 0.5, 2, and 4 10 21 cm" 2 ). Yellow contours mark Af ou t 
(Eq.[3]) and N\> (Eq.[K)]). The yellow square marks the column 
density peak and the reference position listed in Table[5] 
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Fig. 26. Dust temperature map of CB 17 (color scale), overlaid 
with contours of the hydrogen column density (white: 10 21 - 
masked by the yellow N out contour - and 10 22 cm" 2 , black: 2, 4, 
6, 8, and 20 1 21 cm" 2 ). Yellow contours mark N out (Eq.[5]) and 
N\> (Eq.[l0]). The yellow square marks the column density peak 
and position of CB 17-SMM. The asterisk marks the position of 
CB 17-IRS (Table[5j see also |Launhardt et al.|2010| ). 
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Fig. 25. Dust temperature map of CB 6 (color scale), overlaid 
with contours of the hydrogen column density (white: 10 21 cm" 2 , 
black: 2, 4, and 6 10 21 cm" 2 ). Yellow contours mark A^ut (Eq.j3J) 
and N\> (Eq.[l0]). The dotted yellow ellipse marks the Ni> approx- 
imation described in Sect. |5.3| The asterisk marks the position of 
the embedded YSO. 



-2 -4 

ARA [arcmin] 

Fig. 27. Dust temperature map of CB 26 (color scale), overlaid 
with contours of the hydrogen column density (white: 10 21 and 
10 22 cm" 2 , black: 0.5, 2, 4, 6, and 8 10 21 cm" 2 ). Yellow contours 



mark A^ut (Eq.[5]) and N\> (Eq. 10 referring to the starless core 
SMM2). Dotted yellow ellipses mark the N\> approximations for 
the two subcores described in Sect. |5.3| The yellow square marks 
the column density peak and position of CB 26-SMM2; the as- 
terisk marks the position of CB 26-SMM1 (Table[5]). 
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Fig. 28. Dust temperature map of CB 27 (color scale), over- 
laid with contours of the hydrogen column density (white: 10 21 
and 10 22 cm -2 , black: 2, 4, 6, and 8 10 21 cm -2 ). Yellow contours 
mark N out (Eq.[3]) and Ny (Eq.[lO]). The yellow square marks the 
column density peak and the reference position listed in Table[5] 
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Fig. 30. Dust temperature map of CB 68 (color scale), overlaid 
with contours of the hydrogen column density (white: 10 21 and 
10 22 cm" 2 , black: 0.5 - masked by yellow N out contour - 2, 4, 
6, 8, and 20 1 21 cm" 2 ). Yellow contours mark N out (Eq.[5]) and 
N\> (Eq.[l0]). The dotted yell ow ellipse marks the N\> approxi- 
mation described in Sect. 53 The asterisk marks the position of 
the protostar and the reference position listed in Table[5] 
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Fig. 29. Dust temperature map of BHR 12 (color scale), over- 
laid with contours of the hydrogen column density (white: 10 21 - 
masked by yellow N out contour - and 10 22 cm -2 , black: 0.5, 2, 4, 
6, 8, and 20 1 21 cm -2 ). Yellow contours mark N out (Eq.j3J) and 
N\> (Eq.[l0]). The dotte d ye llow ellipse marks the N\> approxima- 
tion described in Sect. |5.3| The two asterisks m ark the positions 
of embedded sources SMM1 and SMM2 (see |Launhardt et aT] 
|20T0] ). 
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Fig. 31. Dust temperature map of B 68 (color scale), overlaid 
with contours of the hydrogen column density (white: 10 21 and 
10 22 cm" 2 , black: 0.5, 2, 4, 6, 8 and 20 1 21 cm" 2 ). Yellow con- 
tours mark A^ut (Eq.[5} and Ny (Eq.[l0]). The yellow square 
marks the column density peak and the reference position listed 
in Tabled 
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Fig. 32. Dust temperature map of CB 130 (color scale), overlaid 
with contours of the hydrogen column density (white: 10 21 and 
10 22 cm -2 , black: 0.5, 2 - masked by yellow N out contour - 4, 6, 
8, and 20 10 21 cm" 2 ). Yellow contours mark N out (Eq.[3) and N\> 
(Eq.[l0]). The yellow circle marks the position of the protostel- 
lar core CB 130-SMM1 with the embedded Class VeLLO. The 
yellow square marks the position of the secondary starless core 
CB 130-SMM2. The asterisk marks the position of the Class I- 
II YSO CB 130-IRS, which most likely associated with CB 13 0, 
but not embedded in its dense core (see Launhardt et al.|2010] ). 




-2 
ARA [orcmin] 

Fig. 33. Dust temperature map of B 335 (color scale), overlaid 
with contours of the hydrogen column density (white: 10 21 and 
10 22 cm" 2 , black: 0.5, 2, 4, 6, 8, and 20 10 21 cm" 2 ). Yellow con- 
tours mark N out (Eq.[3]) and N\> (Eq.[K)]). The asterisk marks 
the position of the protostar and the reference position listed in 
Tabled 
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Fig. 34. Dust temperature map of CB 230 (color scale), overlaid 
with contours of the hydrogen column density (white: 10 21 and 
10 22 cm -2 , black: 2, 4, 6, 8, and 20 1 21 cm" 2 ). Yellow contours 
mark N out (Eq.[3]) and N\> (Eq.[l0]). The two asterisks mark the 
positions of the embedded YSOs (see Launha rdt et al.|2010 ). 
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Fig. 35. Dust temperature map of CB 244 (color scale), over- 
laid with contours of the hydrogen column density (white: 10 21 
and 10 22 cm" 2 , black: 0.5, 2, 4, 6, 8, and 20 1 21 cm" 2 ). Yellow 
contours mark A^ ou t (Eq.[3| and N\> (Eq.[l0]) for the starless core 
CB 244-SMM2. The dotted yellow ellipse marks the N v approx- 
imation for CB244-SMM1 described in Sect.[53] The yellow 
square marks the column density peak of CB 244-SMM2 and 
the reference position listed in Table[5] The asterisk marks the 
position of of the protostar CB244-SMM1. 
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Appendix A: Images on large-scale morphology 

These figures are discussed in Sect.|5.8| 
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Fig.A.l. CB4: Visual (red) DSS2 image with SPIRE 500 //m 
continuum contours overlaid. The red arc marks the warm rim 
ofCB4. 
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Fig. A.3. CB 17: Visual (red) DSS2 image with SPIRE 500 //m 
continuum contours overlaid. The red arc marks the warm rim 
ofCB17. 
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Fig. A.4. CB 26: Visual (red) DSS2 image with SPIRE 500 yum 
continuum contours overlaid. The red arc marks the warm rim 
of CB 26. 
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Fig.A.2. CB 6: Visual (red) DSS2 image with SPIRE 500 yum 
continuum contours overlaid. The red arc marks the warm rim 
ofCB 6. 
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Fig. A.5. CB 27: Visual (red) DSS1 image with SPIRE 500 yum 
continuum contours overlaid. The red arc marks the warm rim 
of CB 27. 
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Fig. A.7. CB 68: Visual (red) DSS2 image with SPIRE 500 /mi 
continuum contours overlaid. The red arc marks the warm rim 
ofCB 68. 
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Fig. A.6. BHR 12: Visual (red) DSS2 image with SPIRE 500 yum Fig. A.8. B 68: Visual (red) DSS1 image with SPIRE 500 /mi 
continuum contours overlaid. The tail side is warmer than the continuum contours overlaid. The red arc marks the warm rim of 
head side, but there is nor clear warm rim side. B 68. The B2 IV star Oph is located at a distance from the Sun 

of 134 pc, i.e., its physical distance from B 68 is only « 2.8 pc. 
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Fig. A.9. CB 130: Visual (red) DSS2 image with SPIRE 500 yum 
continuum contours overlaid. The red arc marks the warm rim 
ofCB130. 
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Fig. A.10. B 335: Visual (red) DSS2 image with SPIRE 500 yum 
continuum contours overlaid. The red arc marks the warm rim 
of B 335. 
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Fig.A.ll. CB230: Visual (red) DSS2 image with 
SPIRE 500 yum continuum contours overlaid. The red arc 
marks the warm rim of CB 230. 
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Fig.A.12. CB244: Visual (red) DSS2 image with 
SPIRE 500 yum continuum contours overlaid. The red arc 
marks the warm rim of CB 244. 



